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I.  INTRODUCTION 


It  has  been  shown  by  Moses, ^  Klarman  et  al.,(2)  Spadaccini  and  Pelmas,^) 
Kinney  and  Lombard,^)  and  others  that  adding  water  to  fuel  in  the  form  of  emulsions 
can  reduce  soot  in  gas-turbine  combustors.  Callahan  et  al.(5)  have  also  demonstrated 
soot  reductions  in  diesel  engines  using  emulsified  fuels.  Several  reasons  have  been 
suggested  for  the  soot  reduction  including:  (a)  a  reduction  in  gas-phase  pyrolysis 
reactions  due  to  the  lowered  combustion  temperature  resulting  from  the  increased 
heat  used  to  vaporize  water  or  an  increase  in  endothermic  reactions,  (b)  a  reduction  in 
liquid-phase  pyrolysis  due  to  lowered  drop  temperatures,  (c)  an  increase  in  soot 
precursor  oxidation  due  to  an  increased  OH  concentration,  or  (d)  an  improvement  in 
fuel/air  mixing  resulting  from  "microexplosions"  caused  by  superheating  of  the  water 
in  the  fuel  drop.  Which  mechanism  dominates  is  still  in  question. 

Recently  published  papers  by  Gollahalli  et  al/6)  and  Rao  and  Bardon^)  show 
opposite  effects  on  flame  temperature  resulting  from  emulsification  for  the  two 
different  flames  studied,  and  attribute  soot  reduction  to  items  (d)  and  (a)  above, 
respectively.  Although  soot  reductions  for  emulsions  observed  in  full-scale  tests  have 
often  been  attributed  to  microexplosions,  it  has  never  been  possible  to  observe  them  in 
these  tests.  Microexplosions  can  occur  when  water  is  contained  in  a  fuel  drop  which 
has  a  higher  boiling  point  than  water.  During  the  heatup  of  such  a  drop,  the  water  may 
either  be  brought  to  the  surface  to  evaporate,  or,  if  that  transport  process  is  too  slow 
and  the  drop  exceeds  the  boiling  point  of  water,  the  superheated  water  may  violently 
boil  and  shatter  the  drop  into  many  smaller  drops.  These  smaller  drops  enhance  mixing 
and  reduce  the  probability  of  soot  formation.  Disruptive  burning  has  been  observed  in 
single-drop  laboratory  experiments^^,  10),  but  the  fuels  and/or  air  conditions  have 
been  chosen  to  maximize  the  probability  of  microexplosions  and  are  not  representative 
of  the  environment  in  a  gas-turbine  combustor.  For  example,  the  high  relative 
drop/gas  velocities  in  a  gas-turbine  combustor  promote  internal  circulation  in  the  fuel 
drops  which  reduce  the  probability  of  disruptive  burning  by  transporting  the  more 
volatile  water  to  the  surface.^)  in  one  single-drop  experiment,  Lasheras  et  al/1Q) 
limited  internal  circulation  by  keeping  the  Reynolds  number  of  the  drops  to  less  than 
0.5,  while  in  a  gas  turbine  the  same  size  drop  would  have  an  initial  Reynolds  number 
about  two  or  three  orders  of  magnitude  larger. 


One  of  the  goals  of  the  work  reported  here  was  to  examine  the  evaporation  of 
emulsified  fuels  sprayed  from  a  pressure  jet  atomizer  into  turbulent  high-temperature 
high-pressure  air,  and  look  for  evidence  of  microexplosions  of  emulsified  drops. 
Comparative  evaporation  rates  of  emulsified  and  neat  fuels  were  determined  using  a 
two-phase  probe,  and  drop-size  distributions  were  determined  from  the  forward  angle 
light  diffraction  pattern  using  a  Malvern  Model  2200  modified  for  use  in  evaporating 
sprays.  Measurements  were  made  at  air  temperatures  to  700  K  and  pressures  to  586 
kPa  using  both  Jet-A  and  a  single-component  fuel,  both  in  err  if ied  and  neat  form. 
Modifications  required  for  drop-size  measurements  in  flames  <-  also  described. 

In  addition  to  the  interest  in  microexplosions,  another  t  of  this  program  has 
been  to  develop  a  better  understanding  of  fuel  spray  evaporat.  oy  coupling  the  drop- 
size  and  fuel  vapor  concentration  measurements  with  a  computer  model  of  spray 
evaporation.  As  indicated  in  Faeth's^^  review  article,  there  are  a  number  of  spray 
models  which  have  been  written,  but  the  verification  data  is  very  sparse.  Data 
generated  in  this  program  directly  address  this  void. 

This  report  describes  the  high-temperature,  high-pressure  spray  test  facility,  the 
operation  and  modifications  to  the  Malvern  drop-sizing  instrument,  the  two-phase 
probe,  and  the  computer  spray  model  that  were  employed  in  this  research  program. 
The  results  obtained  using  these  diagnostics  to  study  fuel-spray  evaporation  for 
emulsified  and  non-emulsified  fuels  are  presented.  Five  publications/presentations 
resulted  from  this  work  and  are  referenced  in  Appendix  B.  These  papers  include 
descriptions  of:  (1)  a  calibration  procedure  for  the  Malvern  instrument  which 
significantly  improves  accuracy  and  is  now  becoming  widely  accepted,  (2)  several 
modifications  to  the  Malvern  which  allow  it  to  be  used  in  evaporating  sprays,  and  (3) 
several  papers  which  show  results  from  these  measurements  and  model  predictions. 


n.  EXPERIMENTAL  APPARATUS 


Four  principal  pieces  of  experimental  apparatus  have  been  assembled  or  pur¬ 
chased  for  this  program.  They  are:  (1)  the  high-temperature  high-pressure  test 
section  and  ^isc-in-duct  combustor,  along  with  the  air  factory,  nozzles,  and  other 
associated  hardware,  (2)  the  Malvern  laser-diffraction,  drop-sizing  instrument  as 
modified  for  these  experiments,  (3)  the  two-phase,  or  phase-discriminating  probe  used 
for  fuel  vapor  measurements,  and  (4)  the  computer  model  for  spray  evaporation. 

2.1  High-Temperature  High-Pressure  Test  Section 

Almost  all  of  the  experiments  described  in  this  report  were  spray  evaporation 
studies  without  combustion,  conducted  in  uniform-temperature  fully-developed  turbu¬ 
lent  flow  in  a  constant  area  duct.  A  few  experiments  were  attempted  in  burning  fuel 
sprays,  but  the  data  were  of  poor  quality  and  difficult  to  interpret.  The  same  test 
section  was  used  for  the  non-combusting  and  combusting  experiments  except  that  a 
disc  coaxial  with  the  circular  duct  was  inserted  for  he  combustion  experiments  to 
provide  a  recirculation  zone  for  flame  stabilization. 

Unvitiated  compressed  air  is  supplied  to  the  test  ceil  by  three  compressors  and  a 
gas-fired  preheater  which  provide  up  to  1.1  kg/sec  at  pressures  ranging  from  138  kPa 
(1.36  atm)  to  1620  kPa  (16.0  atm)  and  temperatures  up  to  1090  K. 

The  test  section  shown  in  Figure  1,  consists  of  a  type  316  stainless  steel  pipe, 
16.83  cm  OD,  12.50  ID,  2.17  cm  wail  thickness,  with  quartz  windows  on  each  side 
having  a  clear  aperture  of  9.2  cm  by  6.7  cm.  A  smooth  transition  section  5.5  cm  long 
is  provided  on  each  end  of  the  windows,  as  shown  in  Figure  2,  to  reduce  recirculation 
zones.  Each  window  is  purged  by  air  on  both  top  and  bottom.  Air  is  supplied  to  each 
window  by  a  total  of  10  tubes  of  nominal  0.635  cm  OD  and  then  diffused  and  exhausted 
by  a  convergent  section  parallel  to  the  window  surface. 

For  the  evaporation  tests,  a  fuel  nozzle  is  mounted  about  even  with  the  upstream 
edge  of  the  windows.  The  nozzle  is  supported  with  a  12.7  mm  OD  stainless  steel  tube 
which  follows  the  test  section  centerline  for  a  distance  of  about  25  cm  upstream  from 
the  nozzle  and  then  makes  a  right  angle  bend  and  exits  from  the  test  section.  The 
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12.7  mm  OD  tube  carries  cooh,.g  water  and  a  separate  parallel  6.35  mm  OD  stainless 
steel  tube  is  the  water  return  line.  Another  6.35  mm  OD  stainless  steel  tube  is 
mounted  inside  the  12.7  mm  water  cooling  tube  and  carries  fuel  to  the  nozzle.  Fuel 
temperature  is  regulated  by  controlling  the  water  temperature. 

For  the  spray  tests  with  combustion,  a  7.22  cm  diameter  disk  holding  the  fuel 
nozzle  in  the  center  and  an  igniter  is  mounted  in  the  center  of  the  duct  as  shown  in 
Figure  3.  The  disk  stabilizer  follows  the  design  by  Mellor  and  coworkers  at  Purdue, 
and  occupies  1/3  of  the  area  of  the  duct.  Control  of  the  fuel  temperature  is  the  same 
as  that  described  above  for  the  evaporation  tests.  For  some  experiments,  the  flame 
was  also  stabilized  without  the  disc  by  using  two  rectangular  flame  holders,  6.4  mm 
thick,  25.4  mm  high,  and  50mm  long,  placed  25  mm  apart  at  the  upper  and  lower  edges 
of  the  spray,  42  mm  downstream  of  the  nozzle. 

The  two  nozzles  used  are  simplex  pressure  atomizers  producing  a  hollow  cone 
spray.  The  model  type  was  Delavan  Corp.  WDA  45°  hollow  cone  with  nominal  flow 
rates  of  3.79  liters/hour  (1  gallon/hour)  and  11.3  liters/hour  (3  gallons/hour)  at  862  kPa 
(125  psi)  differential  pressure  using  water. 

The  operation  of  the  test  cell  is  monitored  with  a  Hewlett-Packard  9820 
programmable  calculator  which  is  coupled  to  a  50  channel  scanner.  A  test  report 
including  the  average  and  statistical  variation  of  ail  the  important  parameters  is 
available  immediately  after  a  test.  The  sensing  systems  consist  of  strain-gage 
pressure  transducers,  thermocouples,  and  turbine  flowmeters. 

2.2  Malvern  Drop-Sizine  Instrument 


Drop  size  data  were  obtained  with  a  Malvern  Model  2200  Particle  Sizer  based  on 
principles  of  forward-light  scattering  or  Fraunhofer  diffraction.  When  illuminated  by 
a  beam  of  monochromatic,  coherent,  collimated  light  from  a  HeNe  laser,  the  smaller 
drops  diffract  light  at  larger  angles  to  the  optical  axis  than  the  larger  drops  and  a 
unique  diffraction  pattern  is  formed.  Detection  is  accomplished  with  a  30  annular  ring 
set  of  solid  state  detectors.  Detector  outputs  are  multiplexed  and  the  data  signal- 
averaged  with  a  Commodore  PET  computer.  A  computer  routine  proprietary  to 
Malvern  is  used  to  interpret  the  light  scattering  pattern  of  the  polydisperse  drop 


systems.  The  drop-size  data  are  available  as  either  a  set  of  two  parameters  defining  a 
Rosin-Rammler,  log-normal,  or  normal  distribution,  or  as  a  histogram  of  15  size 
classes  of  drops  without  any  assumption  about  the  shape  of  the  distribution.  Although 
all  of  the  distributions  have  been  tried,  the  Rosin-Rammler  has  generally  provided  the 
best  fit  and  was  used  throughout  this  program.  A  300  mm  focal  length  f/7.3  lens  was 
used  to  collect  the  scattered  light.  The  laser  beam  diameter  was  9  mm  with  a 
Gaussian  intensity  distribution  truncated  at  the  edge  by  the  9  mm  aperture. 

The  standard  Malvern  system  is  not  suitable  for  use  with  evaporating  sprays  due 
to  beam-steering  problems;  further  complications  are  introduced  if  the  spray  is 
burning  because  the  laser  signal  must  be  separated  from  the  flame  radiation. 
Modifications  were  made  to  the  Malvern  system  to  address  both  of  these  problems. 

The  laser  beam-steering  problems  at  elevated  air  temperatures  were  caused  by 
refractive  index  gradients  due  to  both  temperature  gradients  and  fuel  vapor  concen¬ 
tration  gradients.  The  problem  is  identical  to  "heat  waves"  observed  over  buildings  in 
the  summer.  The  result  is  that  the  unscattered  portion  of  the  laser  beam  of  the  drop 
sizer,  which  normally  is  focused  onto  the  central  detector,  is  caused  to  wander  onto 
some  of  the  other  nearby  detectors  which  are  used  to  measure  light  scattered  by  the 
drops  in  the  spray.  Since  large  drops  scatter  light  at  the  smallest  angle  and  their 
scattered  light  is  focused  onto  the  inner  detectors,  the  extra  signal  intensity  due  to 
the  beam  steering  is  interpreted  by  the  instrument  as  more  large  drops  than  actually 
exist.  Fortunately,  the  additional  signals  due  to  beam  steering  are  typically  limited  to 
the  inner  six  or  eight  of  the  thirty  detectors  in  non-burning  evaporating  sprays  early  in 
the  evaporation  process.  The  detector  radii  increase  rapidly  in  going  away  from  the 
center  spot,  so  that  the  beam  steering  perturbations,  which  affect  the  first  couple  of 
detectors  at  any  elevated  temperature  must  become  much  larger  in  magnitude  to 
progress  outward  through  the  larger  radii  detectors. 

The  signals  from  the  30  detectors  are  paired  together  for  adjacent  detectors  to 
form  15  data  which  are  used  to  determine  the  two  parameters  specifying  the  Rosin- 
Rammler  distributions.  Mathematically,  this  represents  an  overdetermined  set  of  15 
data  to  determine  2  coefficients,  and  some  of  these  data  are  unnecessary  except  to 
improve  the  statistics  in  determining  the  fit.  Thus,  in  some  cases  the  data  from  the 
inner  detectors  (up  to  8)  may  be  ignored  and  a  reasonable  fit  to  the  Rosin-Rammler 


distribution  obtained.  This  approach  has  been  used  successfully  and  is  described  in 
more  detail  by  Dodge  and  Cerwin.^2)  This  modification  allows  the  use  of  this  type  of 
automated  sizing  instrumentation  in  elevated  temperature  and  pressure  sprays  for  the 
first  time. 

Additional  changes  were  required  for  burning  sprays.  Three  modifications  were 
made  to  discriminate  against  flame  radiation.  First  the  laser  beam  was  chopped  at 
667  Hz  and  the  differential  signal  at  each  detector  ring  was  sampled  for  20  cycles. 
This  necessitated  slowing  the  multiplexing  rate  in  the  electronics  from  about  1.3 
ms/channel  to  30  ms/channel.  It  was  also  necessary  to  sample  on  the  "clean"  part  of 
the  signal  when  the  beam  was  completely  blocked  or  completely  transmitting  as  the 
chopper  blade  edge  produces  a  large  spurious  signal  when  passing  through  the  laser 
beam  due  to  its  own  diffraction  pattern.  The  second  modification  was  the  use  of  a 
3  nm  band  pass  interference  filter  centered  at  the  HeNe  laser  wavelength  of  632.8  nm 
in  front  of  the  lens  on  the  drop  sizer.  All  measurements  were  made  with  the  300  mm 
focal  length  lens  which  has  a  maximum  effective  acceptance  angle  of  about  2.7°  from 
the  normal  axis  for  light  reaching  the  detector.  The  transmission  of  the  filter  changes 
less  than  one  percent  for  rays  at  this  angle  or  less  relative  to  the  normal  axis  through 
the  filter.  However,  if  the  100  mm  or  63  mm  focal  length  lenses  (the  other  standard 
sizes)  were  used,  the  acceptance  angle  increases  to  about  8.1°  or  12.8°,  respectively, 
and  wider  band  pass  filters  would  be  necessary.  The  third  modification  involves  the 
use  of  an  aperture  of  diameter  32  mm  located  about  130  mm  from  the  lens  to  block 
radiation  from  the  flame  which  is  not  within  the  scattering  volume  of  the  laser  beam. 
These  modifications  are  discussed  in  more  detail  by  Dodge  and  Cerwin.^2) 

The  modifications  to  the  instrument  to  discriminate  against  flame  radiation  were 
successful  in  allowing  the  light  scattered  by  the  drops  to  be  detected  while  discrimi¬ 
nating  against  flame  radiation  which  was  orders  of  magnitude  higher  in  intensity. 
However,  for  the  conditions  examined,  the  extreme  temperature  gradients  caused  such 
large  beam  steering  excursions  that  reasonable  drop-size  data  were  not  obtained  for 
burning  sprays. 

In  addition  to  the  various  modifications  already  discussed,  a  calibration  proce¬ 
dure  has  been  developed  and  described  in  Reference  3  of  Appendix  B  which  substan¬ 
tially  improves  the  accuracy  of  the  Malvern  instrument.  This  procedure  involves 
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determining  the  responsivity  of  each  of  the  30  detectors  using  a  uniform  light  source. 
After  using  this  procedure,  the  instrument  response  agrees  with  calibration  standards 
within  the  accuracy  of  those  standards. 

Z3  Two-Phase  Probe 

The  temperature-controlled  two-phase  probe  is  shown  in  Figure  4.  The  operation 
of  various  two-phase  probe  designs  has  been  described  by  Collins/13)  McVey  et  al./14>) 
and  Wadleigh  and  Oman/13)  The  concept  is  to  transport  gases  past  a  sample  probe 
(the  vapor-only  probe)  at  a  high  velocity  and  to  pull  a  small  amount  of  vapor  and  air 
out  of  the  high  speed  stream  for  analysis.  The  streamlines  into  the  vapor-only  tube 
bend  so  sharply  that  drops  larger  than  a  few  micrometers  cannot  make  the  sharp  turn, 
and  centrifugal  force  causes  them  to  be  separated  from  the  vapor  and  air  sample.  This 
sample  can  then  be  analyzed  for  fuel  vapor  concentration.  A  second  tube  concentric 
with  the  vapor-only  tube,  but  larger,  is  shown  in  Figure  4.  It  is  called  the  boundary- 
layer  suction  tube,  and  it  prevents  drops  from  impinging  on  the  vapor-only  tube  and 
then  being  sucked  in.  Guidelines  for  operation  of  the  boundary-layer  suction  tube  are 
given  by  Wadleigh  and  Oman/13)  Flow  was  established  in  both  the  boundary  layer 
suction  tube  and  the  main  probe  without  pumps,  but  rather  using  the  pressure 
differential  between  the  test  section  (448  kPa)  and  atmospheric  pressure.  A  valve  was 
used  in  the  boundary  layer  tube  to  restrict  flow  and  was  adjusted  to  give  a  minimum 
reading  (actually  zero)  for  the  vapor-only  tube  gases  while  sampling  a  hexadecane 
spray  at  room  temperature.  Flow  through  the  main  probe  was  not  restricted  to 
achieve  a  maximum  velocity  at  the  vapor/drop  separation  point.  The  main  tube  is  4.93 
mm  id  and  610  mm  long.  Computations  indicated  that  the  flow  should  be  friction- 
choked,  or  close  to  it,  in  the  main  probe,  with  an  entrance  mach  number  of  about  0.40. 
Differential  measurements  of  the  static  and  stagnation  pressure  indicated  an  entrance 
mach  number  of  about  0.33. 

The  boundary  layer  suction  tube  is  2.57  mm  id.  The  vapor-only  tube  is  located 
25  mm  from  the  entrance  to  the  main  probe  and  is  1.08  mm  id.  The  probe  is 
maintained  at  about  478  K  with  recirculating  silicon  oil  which  is  temperature- 
controlled  by  electric  heaters  and  water  cooling.  The  two-meter  long  sample  line  was 
teflon  and  heated  to  about  505  K.  The  analysis  of  the  vapor-only  sample  was  made 
with  a  heated  (422  K)  flame  ionization  detector. 
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A  computer  model  based  on  an  array  of  different-sized,  noninteracting  drops  has 
been  developed  to  study  the  heat-up,  final  "wet-bulb"  drop  temperature,  evaporation 
rates,  trajectories  and  resulting  size-distribution  parameters  at  locations  downstream 
from  the  nozzle.  The  trajectory  part  of  the  model  should  be  considered  as  a 
straightforward  and  simple  one  based  on  an  average  drop  path  (neglecting  turbulent 
statistical  fluctuations)  in  a  one-dimensional  flow  field.  This  corresponds  to  the 
experimental  situation  being  studied.  The  thermodynamic  part  of  the  model,  however, 
is  very  detailed  and  computes  the  transient  heat-up  and  final  drop  temperature  in 
greater  detail  than  many  spray  models.  The  environment  surrounding  the  drop  is 
calculated  based  on  the  properties  of  both  the  fuel  vapor  and  air  rather  than  air  alone 
as  in  some  spray  models.  The  convective  effects  on  evaporation  are  included  and 
constantly  updated  as  the  different  size  drops  are  decelerated  to  the  air-stream 
velocity.  The  effects  of  temperature  and  pressure  on  air  and  fuel  properties  is 
included,  and  the  steady-state  "wet-bulb"  temperature  is  iteratively  calculated  using 
the  equations  for  the  vapor  pressure  and  temperature-dependent  fuel  and  air 
properties.  This  model  assumes  uniform  drop  temperatures  (infinite  heat  transfer 
rates  within  the  drop)  and  single-component  fuels. 

Some  spray  models  have  predicted  the  change  in  drop-size  distribution  for  a 
quiescent  spray  as  a  function  of  time,  but  results  from  this  model  show  the 
considerable  distinction  between  drop-size  distributions  as  measured  at  a  given 
location  along  the  nozzle  axis,  at  which  drops  arrive  after  different  transit  times 
(depending  on  diameter),  as  contrasted  with  the  variation  with  time.  For  example,  in 
this  experiment,  the  spray  initially  has  a  higher  axial  velocity  than  the  air  and  the 
small  drops  decelerate  more  rapidly  than  the  large  ones,  resulting  in  a  higher 
concentration  in  the  measurement  volume  for  the  small  drops  and  a  longer  time-of- 
arrival.  These  effects  must  be  included  in  any  model  which  is  intended  to  be  used  for 
comparison  with  experimentally  measured  changes  in  size  distribution. 

In  developing  a  model  which  can  be  used  to  compare  with  experimental  results,  it 
is  necessary  to  consider  the  interaction  of  the  measuring  device  with  the  spray.  The 
added  weighting  of  the  small  drops  due  to  their  more  rapid  deceleration  has  been 
discussed.  In  addition,  the  small  drops  are  caught  up  in  the  air  stream  and  tend  to 
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populate  the  center  of  the  spray  while  the  large  drops  maintain  their  initial  velocity 
longer  and  tend  to  fill  the  outer  edges  of  the  spray.  This  disparity  is  magnified  by  the 
fact  that  at  the  exit  of  a  hollow-cone  swirl  nozzle  (the  type  used  in  these  tests)  the 
majority  of  the  fluid  consists  of  larger  drops  in  the  outer  cone  of  the  spray,  while 
smaller  drops  fill  in  the  less  dense  spray  in  the  center  of  the  cone.  Thus,  it  is 
necessary  to  model  these  inhomogeneities  in  size  distribution  through  the  spray  and 
also  the  location,  diameter,  and  intensity  distribution  of  the  laser  beam  relative  to  the 
spray.  At  this  time,  the  model  does  account  for  the  different  trajectories  of  the  large 
and  small  drops,  but  assumes  that  all  the  spray  is  initially  at  a  fixed  cone  angle.  The 
drop-size  distribution  is  calculated  for  a  complete  cross  section  of  the  spray  at  a  given 
axial  location,  without  regard  for  the  laser  beam  sampling  volume.  The  model  also 
assumes  that  spray  breakup  is  complete  and  no  secondary  atomization  or  drop 
agglomeration  occurs  beyond  some  point  from  the  nozzle  exit.  An  initial  size 
distribution  at  that  point  is  an  input  to  the  model.  Thus,  the  model  is  not  applicable 
very  close  to  the  nozzle  (  &10mm). 

Results  from  the  model  are  presented  in  the  following  sections,  while  the 
detailed  mathematics  of  the  model  are  presented  in  Appendix  A.  A  general  flow  chart 
of  the  model  is  presented  in  Figure  Al.  Most  of  the  model  has  been  constructed  from 
pieces  of  existing  models.  The  thermodynamic  part  of  the  model  follows  closely  the 
work  of  Chin  and  Lefebvre  at  Purdue  University,  while  the  aerodynamic  equations  for 
drag,  drop  velocity,  etc.  are  based  on  a  model  developed  at  the  University  of  Sheffield 
by  Swithenbank,  Boyson,  Ayers  and  others.  (See  Appendix  A  for  detailed  references). 

The  capabilities  of  the  model  can  probably  be  best  understood  by  examining  some 
typical  output  data  and  also  comparing  the  model  results  with  experimental  results. 
Some  information  available  from  the  model  is  shown  in  Table  2  which  represents  a 
spray  of  n-hexadecane  into  air  at  450  K  and  4.42  atm.  The  motion  of  the  fuel  spray 
and  air  are  illustrated  in  Figure  5.  The  air  velocity  was  assumed  to  be  a  uniform  one¬ 
dimensional  axial  flow  of  8.54  m/s,  while  the  initial  fuel  velocity  was  20.76  m/s  at  an 
angle  of  22.5°  to  the  axis  (45°  cone  angle)  and  a  fuel  temperature  of  311  K.  The 
computed  drop  size  data  assuming  a  Rosin-Rammler  distribution  are  shown  on  the  top 
line  and  include  the  two  Rosin-Rammler  parameters,  x  and  N  which  specify  a  size  and 
a  width  for  the  distribution,  and  the  Sauter  mean  diameter  (SMD)  which  is  an 
"average"  drop  size.  By  definition,  if  all  drops  in  the  spray  were  equal  in  size  to  the 
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SMD,  than  the  volume-to-surface-area  ratio  of  the  imaginary  monodisperse  spray 
would  be  equal  to  the  actual  spray.  This  ratio  is  chosen  because  surface  area  regulates 
evaporation  rates  while  the  volume  indicates  the  quantity  of  fuel.  The  initial  drop 
sizes  are  shown  in  column  1,  and  the  initial  relative  population  of  drops  in  each  class 
were  computed  to  correspond  exactly  to  a  Rosin-Rammler  distribution  specified  as  an 
input  to  the  program.  (Those  values  were  x  =  52  and  N  =  1.75  for  the  example  shown). 
Because  of  the  evaporation  and  the  more  rapid  deceleration  of  the  small  drops 
(resulting  in  an  increased  weighting  for  the  small  drops),  the  size  distribution  changes 
at  positions  downstream  from  the  nozzle.  The  drop-size  distribution  no  longer 
corresponds  exactly  to  the  Rosin-Rammler  distribution,  and  the  degree  to  which  it 
does  fit  that  distribution  is  given  by  the  correlation  coefficient.  For  the  cases 
examined,  the  correlation  coefficient  has  been  in  excess  of  0.97,  but  under  some 
conditions  it  could  be  worse.  These  distribution  parameters  are  calculated  from 
information  shown  for  individual  drop  size  classes  in  the  lower  part  of  Table  2. 

The  initial  drop  sizes  (in  Mm)  are  shown  in  Column  1,  while  the  instantaneous  drop 
sizes  are  shown  in  column  2  at  a  distance  from  the  nozzle  of  16  mm  as  shown  in 
column  4  and  a  time  (in  seconds)  as  specified  in  column  3.  When  the  drops  evaporate 
completely,  the  time  and  location  are  frozen  at  the  point  where  the  drop  disappeared 
and  those  drops  are  not  used  in  calculating  the  Rosin-Rammler  parameters.  Note  that 
in  this  example,  the  model  calculates  spray  data  at  a  given  axial  distance  from  the 
nozzle  (X  DROP)  and  not  at  a  given  time  from  injection.  The  model  can  also  use  time 
as  a  target.  The  smaller  drops  take  longer  to  reach  the  specified  location  due  to  their 
lower  velocity  as  specified  in  column  7  by  U-VEL  in  m/s.  The  smaller  drops  decelerate 
to  the  air  stream  velocity  of  8.54  m/s  more  rapidly  than  the  large  drops.  The 
coordinate  system  is  cylindrical  with  the  axial  coordinate  being  x  with  velocity  u,  the 
corresponding  radial  components  y  and  v,  and  the  corresponding  angular  components 
(assumed  zero  presently)  z  and  w.  Distances  are  in  meters  and  velocities  in 
meters/seconds. 

Note  the  difference  in  trajectories  and  velocities  for  the  different  size  classes. 
All  non-evaporated  drops  are  at  the  same  axial  (X  DROP)  location,  but  the  radial 
location  (Y  DROP)  of  the  small  drops  is  less  than  the  large  ones  due  to  their  being 
caught  up  in  the  air  stream  sooner  as  shown  by  the  axial  (U-VEL)  and  radial  (V-VEL) 


velocities.  The  importance  of  convective  effects  in  heat  and  mass  transfer  is  obvious 
from  these  values,  and  hence,  the  shortcomings  inherent  in  quiescent  spray  models. 

The  instantaneous  drop  temperature  (FTEMP)  is  shown  in  the  next  to  the  last 
column.  The  initial  fuel  temperature  was  311  K  and  the  "wet-bulb"  final  steady-state 
drop  temperature  is  about  433  K  at  these  conditions.  The  normal  boiling  point  of  n- 
hexadecane  is  560  K.  Note  that  the  small  drops  heat  up  and  evaporate  almost 
instantaneously  while  the  "uniformly  mixed"  temperature  of  the  largest  drops  is 
predicted  to  be  close  to  the  initial  fuel  temperature.  Due  to  finite-rate  heat  transfer, 
the  surface  temperature  of  the  large  drops  would  be  higher,  but  experimental  evidence 
is  lacking  to  predict  temperature  profiles  within  the  drop.  Complex  models  have  been 
developed  to  predict  these  profiles,  but  they  are  too  cumbersome  for  this  effort  and 
generally  lack  verification. 

Note  the  importance  of  modeling  the  transient  nature  of  the  drop  before  it 
reaches  the  steady-state  temperature.  Of  the  21  drop  size  classes  shown  in  Table  2, 
only  five  size  classes  have  reached  the  steady-state  temperature.  As  the  pressure  and 
convective  effects  are  increased,  the  length  for  the  transient  period  relative  to  the 
steady-state  period  increases^16).  Thus,  models  that  assume  no  mass  transfer  until  the 
drop  reaches  its  steady-state  temperature  are  not  representative  of  many  realistic 
conditions  in  gas  turbines  and  are  insufficient  to  compare  with  experimental  results 
generated  in  this  program. 

The  top  part  of  Table  2  is  computed  from  the  data  for  the  individual  size  classes 
shown  in  the  bottom  of  Table  2.  The  first  line  of  the  middle  section  indicates  the 
fraction  of  the  fuel  which  has  evaporated  at  a  given  axial  location  measured  from  the 
nozzle  exit  (in  meters),  and  the  second  line  indicates  the  fuel  remaining  in  liquid  form. 
The  fuel  vaporization  absorbs  heat  from  the  surrounding  air  in  the  spray  cone  and  the 
resulting  air  temperature  is  shown  on  the  third  line.  Momentum  transfer  from  the 
spray  to  the  air  is  not  included  in  this  model,  but  the  effect  of  lowered  temperature  on 
increasing  air  density  and  reducing  air  velocity  in  the  spray  cone  is  included,  as  shown 
on  the  fourth  line.  The  quantity  of  vaporized  fuel  is  calculated  and  assumed  to  be 
uniformly  distributed  at  a  given  axial  location  over  the  region  shown  in  Figure  5.  The 
fuel  vapor  concentration  expressed  as  a  mass  fraction  is  shown  on  the  fifth  line  and  as 
mole  fraction  on  the  sixth  line.  In  order  to  compare  these  values  with  those  measured 
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by  the  probe,  the  mole  fraction  integrated  over  a  cross-section  of  the  spray  cone  at  a 
given  axial  location  is  given  in  line  seven. 

In  addition  to  the  drop  array  data  just  discussed,  the  model  has  been  used  to 
examine  the  steady-state  "wet-bulb"  temperature  and  the  quiescent,  evaporation  rate 
of  the  fuel,  which  are  both  independent  of  drop  size,  as  a  function  of  air  pressure  and 
temperature.  At  atmospheric  pressure,  the  steady-state  drop  temperature  is  lower 
than  the  normal  boiling  point  due  to  the  cooling  effect  of  energy  used  in  vaporizing  the 
fuel.  As  the  pressure  increases,  the  boiling  point  increases  and  the  steady-state 
temperature  increases.  Obviously  the  fuel  and  air  properties  are  important  in 
determining  this  temperature.  Steady-state  temperatures  for  three  normal  paraffins  - 
octane  (CgHig),  dodecane  (C  12^26^  anc*  hexadecane  (C16H34)  -  were  calculated  over 
a  range  of  conditions  with  the  results  as  shown  in  Figures  6  to  8.  At  low  air 
temperature:,  the  drop  temperatures  converge  to  the  air  temperature.  As  the  air 
temperature  increases,  the  increased  vaporization  cools  the  drop  below  the  air 
temperature,  but  ncreased  pressure  markedly  increases  the  allowable  drop  tempera¬ 
ture.  The  spontaneous  nucleation  temperature  for  water  indicated  in  the  figures  is 
discussed  in  a  later  section.  The  critical  pressure  for  hexadecane  is  1421  kPa  (14.02 
atm)  and  for  dodecane  1820  kPa  (17.96  atm),  and  data  are  not  shown  for  pressures 
higher  than  the  critical  pressure. 

Droplet  evaporation  follows  a  D-squared  law  which  relates  the  drop  diameter  D 
to  the  initial  size  D0  and  the  time  t  by^2\ 

D2  =  D02  -  \t  (1) 

where  the  proportionality  constant  X  is  the  evaporation  constant.  The  mass  transfer 
rate  of  fuel  from  the  drop,  mp,  at  quiescent  conditions  is  related  to  the  evaporation 
constant  by, 


rhp  =  (tt/4)  X Pp  D 


(2) 


where  pf  is  the  fuel  density.  The  evaporation  constant  is  increased  by  convection, 
with  a  correction  factor  given  by  Frossling(^)  for  cases  where  heat  transfer  rates  are 

controlling  of, 


SHE 


xConv.  =  ^Quies.  ^  +  0.276  Reo^  Prg®*^^  ) 


where  Ref)  is  the  Reynolds  number  of  the  gas  relative  to  the  drop  and  Prg  is  the 
Prandtl  number  of  the  gas. 

In  summary,  the  computer  model,  while  not  completely  comparable  with 
experimental  results  in  its  present  form,  does  provide  insight  into  much  of  the 
thermodynamics  affecting  drop  evaporation.  It  also  shows  quantitatively  how  the 
smellier  drops  are  caught  up  in  the  air  stream  much  more  quickly  than  the  larger  ones. 
It  also  indicates  the  importance  of  modeling  evaporation  during  the  transient  heat-up 
phase,  particularly  for  realistic  sprays  at  high  pressures  and  Reynolds  numbers.  It 
shows  the  effect  of  air  pressure  on  the  steady-state  drop  temperature,  and  how 
increasing  air  pressure  leads  to  increased  drop  temperatures  which  enhance  the 
probability  of  initiating  microexplosions. 


in.  EXPERIMENTAL  PROCEDURES  AND  RESULTS 

Experimental  results  include  the  effects  of  fuel  properties,  nozzle  pressure  drop, 
air  temperature,  and  air  pressure  on  the  initial  atomization  of  fuels  and  the  change  in 
drop-size  distribution  during  evaporation.  The  properties  of  the  fuels  used  in  these 
tests  are  given  in  Table  1.  Most  of  these  tests  were  comparative  studies  for  the  fuels 
emulsified  with  20  percent  water  and  the  same  fuels  in  non-emulsified  or  "neat"  form. 
The  20  percent  emulsion  was  chosen  to  accentuate  the  differences  between  the 
emulsified  and  neat  fuel.  The  changes  in  the  experimentally  determined  drop-size 
distributions  with  distance  from  the  nozzle  are  compared  with  values  computed  from 
the  spray  model.  The  spray  model  was  also  used  to  compute  drop  temperatures  to  help 
interpret  the  experimental  results.  Experimental  results  are  also  reported  for  the 
two-phase  probe  measurements  of  the  comparative  vaporization  rates  of  emulsified 
and  neat  fuels. 

3.1  Spray  Drop-Size  Measurements 


3.1.1  General  Characteristics  of  Sprays 

A  comparison  of  the  spray  data  on  centerline  and  at  the  edge  of  the  spray  is 
instructive  in  understanding  the  features  of  fuel  sprays  from  swirl  nozzles.  Typically, 
the  attenuation  of  the  laser  beam  of  the  drop  sizer  by  the  edge  of  the  spray  is  about 
comparable  to  that  for  the  beam  passing  through  the  center  of  the  spray,  although  the 
path  length  through  the  edge  is  much  shorter  than  the  center.  This  can  be  understood 
by  considering  that  in  swirl  atomizers,  the  majority  of  the  fuel  goes  into  larger  drops 
along  the  outer  edge  of  the  cone,  with  am  inner  core  composed  of  smaller  drops  at  a 
much  lower  density.  Figure  9  shows  a  comparison  of  the  characteristic  drop  sizes, 
represented  by  the  Sauter  mean  diameter  (SMD),  measured  along  the  edge  of  the  spray 
and  the  centerline.  These  data  are  for  a  nozzle  similar  in  design  but  larger  in  capacity 
than  for  other  data  in  this  report,  but  the  nozzle  spray  pattern  is  representative. 
These  results  are  consistent  with  the  above  description  of  the  spray  in  that  near  the 
nozzle  the  measurements  are  close  to  the  same  place  (considering  the  9  mm  diameter 
of  the  measuring  beam),  and  SMD's  are  similar,  while  further  downstream  the  edge 
measurements  are  weighted  towards  the  large  drops  and  the  centerline  measurements 
are  more  representative  of  the  overall  distribution.  Further  evidence  for  this  view  of 
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Table  1.  Some  Properties  of  Fuels  Used 


Jet  A  (AL-10112F) 


Viscosity  at  40°C  i  .72  cS 
Specific  Gravity  at  60°F  0.8049 
Heat  of  Combustion  46.904  MJ/kg 
Hydrogen  Content  14.15% 

Flash  Point:  355K 

Boiling  Point  Distribution  (GC) 
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20%  H20/2%  Surfactant/78%  Jet  A  (AL-10112F) 


Viscosity  at  40°C  4.47  cS 
Specific  Gravity  at  60°F  0.8618 


Hexadecane  (95%  purity) 


Viscosity  at  40°C  2.98  cS 
Specific  Gravity  at  60°F  0.7774 
Boiling  Point  (nominal)  560K 


20%  H20/2%  Surfactant/78%  Hexadecane 


DegK 
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Viscosity  at  40°C  5.55  cS 

Specific  Gravity  at  60°F  0.8249 
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FIGURE  9.  COMPARISON  OF  SMD  ON  CENTERLINE  AND  EDGE  OF  SPRAY 


the  spray  structure  is  shown  in  Figure  10  which  compares  the  Rosin-Rammler  N 
parameter  for  the  edge  and  centerline  measurements.  The  N  parameter  is  defined 
by,  (19, 20) 

R  =  exp  (-((d/3N))  (1) 

where  R  is  the  cumulative  volume  fraction  of  drops  larger  than  size  d,  X  is  a  size 
parameter,  and  N  measures  the  width  of  the  distribution.  Large  values  of  N  imply 
narrow  distributions  and  vice  versa.  Thus,  Figure  10  shows  that  close  to  the  nozzle, 
the  measurements  agree,  as  they  should,  while  further  from  the  nozzle  the  edge 
measurements  show  a  much  narrower  distribution  consisting  of  just  the  larger  drops 
(from  Figure  9),  while  the  distribution  through  the  centerline  remains  relatively  broad 
since  both  larger  drops  on  the  edges  and  smaller  drops  in  the  center  of  the  cone  are 
being  sampled.  For  completeness,  it  should  also  be  noted  that  the  SMD  can  be 
computed  from  the  two  Rosin-Rammler  parameters  by(20) 

SMD  =  X/  r(l-l/N)  (2) 

where  T  is  the  gamma  function. 

Further  evidence  of  the  non-homogeneity  of  the  spray  size  distribution  is  shown 
in  Figure  11  for  a  5.68  liter/hr  (1.5  gph)  80°  hollow  cone  nozzle.  This  figure  shows  the 
SMD  and  N  values  at  a  fixed  axial  distance  of  25  mm  at  different  distances  from  the 
centerline.  Again  the  measured  size  distribution  through  the  spray  is  the  smallest  in 
average  size  and  broadest  in  distribution  through  the  center  of  the  spray,  while  the 
edges  are  more  monodisperse  and  much  larger  in  average  size. 

Several  other  comments  can  be  made  concerning  the  trends  in  SMD  and  N  values 
as  a  function  of  distance  from  the  nozzle.  Returning  to  Figure  9,  the  SMD's  on  the 
centerline  decrease  in  going  downstream  from  the  nozzle.  This  is  probably  due  mostly 
to  the  fact  that  close  to  the  nozzle,  the  laser  beam  samples  more  of  the  larger  drops 
along  the  top  and  bottom  edges  of  the  spray,  while  further  downstream  it  is  weighted 
more  heavily  towards  the  smaller  drops  in  the  center  of  the  spray  cone.  Two  other 
effects  contribute  to  this  trend.  First,  when  sampling  very  close  to  the  nozzle, 
( S  12  mm)  the  SMD  is  larger  due  to  the  fact  that  atomization  is  not  complete. 


Second,  in  the  first  few  centimeters  downstream  of  the  nozzle  the  small  drops  slow 
down  to  the  air  stream  velocity  more  rapidly  than  the  large  drops  (see  Table  2), 
increasing  their  relative  population  in  the  laser  beam  sampling  region.  Thus,  in  a  non¬ 
evaporating  spray,  the  SMD  decreases  over  the  first  5  to  10  cm  and  then  levels  out. 

Evaporation  of  the  fuel  drops  has  some  interesting  effects  on  SMD  and  the  N 
parameter.  If  the  spray  were  monodisperse  (single-sized  drops)  then  the  SMD  would 
decrease  with  evaporation.  Sprays  from  pressure  atomizers  are  very  polydisperse 
(broadly  distributed  in  size)  with  N  parameters  on  the  order  of  2.  The  effect  of 
evaporation  is  to  preferentially  evaporate  the  smaller  drops,  in  accordance  with  the 
d^-law,  reducing  the  width  of  the  distribution  and  increasing  the  average  size  for  those 
drops  remaining.  Thus,  at  elevated  temperatures,  it  is  expected  that  in  going 
downstream  from  the  nozzle,  the  SMD  and  N  values  should  both  increase  through  the 
early  stages  of  evaporation.  This  is  predicted  by  the  computer  model  and  has  been 
observed  experimentally,  as  will  be  shown  in  the  following  results. 

Another  way  of  examining  the  fuel  distribution  from  a  spray  nozzle  is  to  measure 
the  quantity  of  fuel  sprayed  at  different  angles  relative  to  the  nozzle  axis.  The  results 
of  such  a  test  for  the  3.79  liter/hr,  45°  hollow-cone  nozzle  are  shown  in  Figure  12  for 
the  cumulative  fraction  of  fuel  within  different  cone  angles  measured  from  the 
centerline.  The  fuel  fraction  at  a  given  angle  is  proportional  to  the  derivative  of  the 
curve,  so  the  inner  10°  on  each  side  of  the  nozzle  axis  (20°  total  included  angle) 
contains  little  fuel,  while  the  majority  falls  between  12°  and  22°  from  the  axis.  The 
nozzle  is  specified  at  45°  total  cone  angle,  corresponding  to  22^°  from  the  axis,  and 
about  86  percent  of  the  fuel  falls  within  the  45°  cone. 

Pressure  drop  across  the  nozzle  has  a  significant  effect  on  the  average  size  and 
width  of  the  spray  distribution  as  shown  in  Figure  13.  Increasing  pressure  drop  causes 
a  reduction  in  SMD  and  a  broadening  of  the  size  distribution.  Also  the  SMD  increases 
slightly  with  viscosity,  the  viscosities  of  the  solvent,  3et-A,  and  the  emulsion  at  40°C 
being  1.00  cS,  1.72  cS,  and  4.47  cS,  respectively. 


FIGURE  12.  ANGULAR  DISTRIBUTION  OF  CUMULATIVE  FUEL  FRACTION 
50  mm  from  3.79  liter /hr  45°  Hollow  Cone  Nozzle 


3.1.2  Effect  of  Elevated  Temperature  and  Pressure  on  Size-Distributions  for 

Emulsified  and  Neat  3et-A 

One  of  the  goals  of  this  program  was  to  study  the  effects  of  emulsification  on 
the  evaporation  of  jet  fuel.  The  effects  of  air  temperature  at  constant  air  pressure  on 
the  fuel  spray  drop-size  distribution  shown  in  Figures  14-18.  Compared  with  non¬ 
evaporating  (relatively)  sprays  at  room  temperature,  it  can  be  seen  that  the  SMD 
increases  during  the  evaporation  process  and  the  higher  the  temperature  the  more 
rapid  the  increase.  The  N  values  also  increase  going  downstream  from  the  nozzle. 
Both  of  these  trends  are  in  agreement  with  the  idea  previously  discussed  that  the 
preferential  evaporation  of  the  smaller  drops  leads  to  a  narrower  distribution  and 
increased  average  size  for  the  remaining  drops.  Both  the  SMD  and  N  values  are 
affected  by  errors  caused  by  beam  steering,  and  these  errors  limit  how  far  downstream 
measurements  were  reported  at  temperatures  of  547  K  and  above.  The  N  values  are 
affected  more  strongly  than  the  SMD's.  These  errors  tend  to  increase  SMD  and 
decrease  N  from  their  actual  values.  The  initial  atomization  of  emulsified  and  neat 
jet-A  is  similar  for  each  of  these  conditions.  Differences  in  initial  drop  sizes  are 
probably  due  mainly  to  differences  in  pressure  drop  across  the  nozzle  and  differences 
in  temperatures  of  the  fuels.  At  450  K,  the  emulsified  and  neat  fuels  show  similar 
trends.  At  589  K  and  644  K,  the  SMD's  for  the  emulsified  and  neat  fuel  diverge  with 
the  emulsified  fuel  producing  smaller  drops  downstream  of  the  nozzle.  All  of  the  SMD 
data  shown  in  Figures  14-18  are  compared  in  Figure  19a,  normalized  to  make 
comparisons  easier.  The  450  K  and  506  K  data  are  normalized  to  the  aveiuge  of  the 
first  two  data  closest  to  the  nozzle.  At  the  higher  temperatures  of  589  K  and  644  K, 
the  evaporation  is  so  rapid  that  only  the  first  points  closest  to  the  nozzle  are  used  to 
normalize  the  data.  The  normalized  SMD's  at  a  fixed  distance  of  16  mm  from  the 
nozzle  are  shown  as  a  function  of  temperature  in  Figure  19b,  with  the  emulsified  fuel 
compared  with  the  neat  fuel.  Figures  19a  and  19b  show  clearly  the  similarity  of  the 
trends  for  drop-size  distributions  during  evaporation  at  the  lower  temperatures  and  the 
divergence  at  the  elevated  temperatures  with  the  emulsified  fuels  producing  smaller 
drops  during  the  evaporation  process.  Similar  trends  were  observed  previously  with  a 
larger  capacity  (11.4  liters/hr)  nozzle/21) 

There  are  at  least  two  possible  explanations  for  the  smaller  drop  sizes  for  the 
emulsified  fuel  relative  to  the  neat  fuel  at  the  higher  temperature  conditions.  One 
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FIGURE  14.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

3ET-A  450  K,  448  kPa 


FIGURE  16.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

3ET-A  547  K,  448  kPa 
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FIGURE  17.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

3ET-A  589  K,  448  kPa 
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FIGURE  18.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

JET- A  644  K,  448  kPa 


FIGURE  19b.  EFFECT  OF  AIR  TEMPERATURE  ON  SMD  OF  EMULSIFIED  AND 
NEAT  JET  A  AT  16  mm.  AIR  PRESSURE  448  kPa 


possible  explanation  for  reduced  drop  sizes  for  emulsified  fuels  relative  to  neat  fuels  is 
microexplosions,  as  previously  discussed.  The  boiling  points  of  many  fuels  are 
sufficiently  high  that  water  could  reach  its  spontaneous  nucleation  point  and  cause  the 
emulsified  drop  to  explode  if  the  water  is  not  transported  to  the  surface  fast  enough 
and  allowed  to  evaporate.  The  spontaneous  nucleation  temperature  of  water  in 
hexadecane  has  been  estimated  by  Avedisian  to  be  about  539  k/^2,23)  ancj  varies  only 
slightly  with  pressure.  This  is  an  upper  limit  for  clean,  degassed  water  in  hexadecane 
to  begin  boiling;  the  actual  boiling  point  for  these  emulsions  could  occur  within  the 
range  of  373  K  to  539  K  depending  on  the  presence  of  bubbles  or  particles  within  the 
water,  and  particularly  at  the  interface  between  the  water  and  the  fuel. 

A  second  possible  explanation  for  the  reduced  average  drop  sizes  of  the 
emulsified  fuel  relative  to  the  neat  fuel  is  that  the  heat  of  vaporization  of  water  is 
about  eight  times  that  of  jet-A,  although  the  boiling  point  of  water  (373  K)  is 
significantly  less  than  that  of  the  "average"  (50  percent)  boiling  point  of  jet-A  (502  K). 
The  higher  heat  of  vaporization  of  water  reduces  evaporation  rates  for  emulsified 
fuels,  but  the  lower  boiling  point  increases  it.  If  the  net  evaporation  rate  of  the 
emulsified  fuels  is  reduced  relative  to  the  neat  fuels,  then  trends  such  as  those  shown 
in  Figures  14-19  could  result. 

Another  similar  set  of  experiments  was  performed  with  emulsified  and  neat  jet- 
A  in  which  the  air  temperature  was  held  constant  and  the  air  pressure  was  varied. 
Law(24)  has  suggested  that  an  increase  in  air  pressure  at  constant  temperature  should 
favor  the  occurrence  of  microexplosions  because  the  steady  state  drop  temperature 
increases  with  air  pressure,  as  was  shown  in  Figures  6-8.  Thus,  for  a  fixed  air 
temperature  of  547  K,  the  spray  drop  sizes  were  determined  as  a  function  of  distance 
from  the  nozzle  at  pressures  of  165  kPa  shown  in  Figure  20,  207  kPa  in  Figure  21, 
345  kPa  in  Figure  22,  and  448  kPa  in  Figure  23.  At  the  lowest  air  pressure,  the 
emulsified  fuel  shows  larger  drops  throughout  the  evaporation  process,  but  at  higher 
pressures,  which  raise  the  steady  state  drop  temperature,  the  SMD's  further  down¬ 
stream  of  the  nozzle  are  smaller  for  the  emulsified  fuels.  These  trends  are  also 
evident  in  Figure  24a  in  which  the  SMD's  for  Figures  20-23  are  all  shown  together. 
Another  way  of  examining  the  pressure  dependence  is  to  compare  the  SMD's  at  a  fixed 
distance  from  the  nozzle  as  a  function  of  pressure.  Such  a  comparison  for  a  distance 
of  26  mm  is  shown  in  Figure  24b,  which  indicates  the  SMD's  after  normalization  to  the 
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FIGURE  20.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

3ET-A,  547  K.  165  kPa 
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N,  Emulsified  3et-A 
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SMD,  Emulsified  3et-A 
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initial  SMD  close  to  the  nozzle.  (The  normalization  is  necessary  because  initial 
atomization  depends  on  air  density.)  The  divergence  of  drop  sizes  at  higher  pressures 
is  consistent  with  the  microexplosion  phenomena. 

In  addition  to  the  effects  of  pressure  on  the  comparative  evaporation  of  neat  and 
emulsified  fuels,  note  that  the  increased  air  density  has  a  significant  effect  on 
improving  the  initial  atomization  close  to  the  nozzle.  This  fact  is  very  significant  in 
that  it  implies  that  nozzle  drop-size  calibrations  performed  at  atmospheric  conditions 
cannot  be  used  to  predict  spray  evaporation  at  elevated  pressure  or  subatmospheric 
conditions  without  correcting  for  this  air  density  effect.  Since  evaporation  times  are 
proportional  to  the  square  of  the  initial  drop  sizes,  the  differences  shown  in  Figure  24a 
for  initial  sizes  represent  even  larger  differences  when  considered  in  terms  of  spray 
evaporation  times. 

3.1.3  Effect  of  Elevated  Temperature  and  Pressure  on  Size-Distribution  For 
Emulsified  and  Neat  Hexadecane 

The  evaporation  of  a  multiple  component  petroleum  fuel  such  as  jet-A  is 
difficult  to  model  even  without  the  addition  of  water.  In  order  to  allow  for  a  more 
quantitative  analysis  of  the  evaporation  process,  experiments  similar  to  those  already 
discussed  for  jet-A  were  conducted  with  emulsified  and  neat  hexadecane.  The 
evaporation  of  these  sprays  was  studied  by  taking  drop-size  measurements  at  a  fixed 
pressure  and  varying  temperature  and  visa  versa.  At  a  fixed  pressure  of  448  kPa  (4.42 
atm),  the  results  are  shown  for  a  series  of  temperatures:  450  K  in  Figure  25,  506  K  in 
Figure  26,  54?  K  in  Figure  27,  589  K  in  Figure  28,  644  K  in  Figure  29,  and  700  K  in 
I  Figure  30.  The  SMD  data  for  Figures  26-30  and  data  for  728  K  are  also  shown  in 

Figure  31a  after  normalization  for  an  overall  comparison  of  behavior.  (The  SMD  data 
for  emulsified  hexadecane  at  506  K  were  identical  to  547  K  and  are  not  shown.)  The 
normalized  SMD's  in  Figure  31a  for  a  distance  of  26  mm  from  the  nozzle  are  replotted 
I  in  Figure  31b  as  a  function  of  air  temperature  to  illustrate  the  temperature 

dependence  more  clearly.  At  450  K  and  506  K,  the  emulsified  and  neat  fuels 
evaporate  in  such  a  way  that  the  drop  sizes  are  very  similar.  At  589  K,  644  K,  700  K, 
and  728  K  the  initial  atomization  is  similar,  but  the  emulsified  fuels  produce  smaller 
I  drops  than  the  neat  fuels  downstream  of  the  nozzle.  These  trends  are  similar  to  those 

observed  for  jet-A  and  are  consistent  with  the  behavior  expected  if  microexplosions 
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^N,  Emulsified  Hexadecane 


FIGURE  25.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

HEXADECANE  450  K,  448  kPa 


Emulsified  Hexadecane 


FIGURE  26.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

HEXADECANE  506  K.  448  kPa 
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FIGURE  27.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

HEXADECANE  547  K,  448  kPa 


♦  N,  Emulsified  Hexadecane 
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FIGURE  28.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

HEXADECANE  589  K,  448  kPa 


o  N,  Neat  Hexadecane 
♦  N,  Emulsified  Hexadecane 


HEXADECANE  644  K,  448  kPa 


Table  3.  Computer  Predicted  Steady-State  Drop 
Surface  Temperature  for  Hexadecane  at  Conditions 
Corresponding  to  Figure  31 
(Pair  =  kPa  (*.42  atm)) 


Air  Temperature  Drop  Temperature 


(K) 

(K) 

506 

469 

547 

488 

589 

503 

644 

520 

700 

535 

728 

541 

Table  4.  Computer  Predicted  Steady-State  Drop 
Surface  Temperature  for  Hexadecane  at  Conditions 
Corresponding  to  Figures  32-35 
(Tair  =  644  K) 


Air  Pressure 

Drop  Temperature 

(kPa) 

(atm) 

(K) 

165 

1.63 

491 

241 

2.38 

502 

345 

3.40 

512 

448 

4.42 

520 

56 


were  occurring.  However,  the  advantage  to  using  a  single-component  fuel  is  that  it  is 
possible  to  predict  the  steady  state  drop  temperature,  and  if  the  assumption  is  made 
that  the  less  volatile  (hexadecane)  component  controls  the  evaporation  rate  for  the 
emulsified  fuel  after  some  initial  period,  then  these  predicted  single-component  values 
should  be  reasonable  estimates  for  the  emulsified  drops.  The  steady-state  drop 
temperatures  for  hexadecane  predicted  by  the  computer  model  described  previously 
are  shown  in  Table  3  for  conditions  corresponding  to  Figure  31a.  The  emulsified  and 
neat  fuels  diverge  in  drop-size  behavior  for  an  air  temperature  of  about  589  K  to  644 
K,  corresponding  to  a  drop  temperature  of  503  K  to  520  K.  The  temperature  for 
microexplosions  to  occur  depends  on  the  purity  and  amount  of  bubbles  at  the 
water/hexadecane  interface,  but  the  upper  limit  for  spontaneous  nucleation  for  a  pure 
bubble-free  system  is  about  539  K  at  atmospheric  pressure,  with  only  a  slight  pressure 
effect/22,23)  js  reasonable  to  expect  that  the  presence  of  bubbles  and  impurities 
may  have  reduced  the  effective  spontaneous  nucleation  temperature  from  a  predicted 
drop  temperature  of  539  K  to  the  range  of  503  K  to  520  K,  in  agreement  with  the 
observed  point  of  divergence  of  emulsified  and  neat  fuels.  The  water/hexa- 
decane/surfactant  system  used  in  these  experiments  was  recirculated  through  an 
ultrasonicator  before  usage  to  maintain  a  well-mixed  blend,  but  this  may  have  added  a 
considerable  number  of  bubbles  into  the  fuel,  lowering  the  effective  spontaneous 
nucleation  temperature  from  the  theoretical  limit. 

As  drop  temperatures  increase,  the  likelihood  of  microexplosions  increases. 
Drop  temperatures  may  be  increased  by  an  increase  in  air  temperature  (at  constant  air 
pressure),  as  discussed  above  and  shown  in  Table  3,  or  by  an  increase  in  air  pressure  (at 
constant  temperature)  as  shown  in  Table  4.  A  set  of  experiments  was  performed  to 
examine  the  effects  of  air  pressure,  and  the  corresponding  steady-state  drop  tempera¬ 
tures,  on  the  comparative  drop-size  distributions  of  emulsified  and  neat  hexadecane  at 
a  constant  air  temperature  of  644  K.  Data  for  four  different  air  pressures  were 
recorded.  The  results  for  165  kPa  (1.63  atm)  are  shown  in  Figure  32,  for  241  kPa  (2.38 
atm)  in  Figure  33,  for  345  kPa  (3.40  atm)  in  Figure  34,  and  for  448  kPa  (4.42  atm)  in 
Figure  35.  The  SMD's  at  a  fixed  distance  of  26  mm  from  the  nozzle  in  Figures  32-35 
are  shown  as  a  function  of  air  pressure  in  Figure  36.  The  SMD's  in  Figure  36  have  been 
normalized  relative  to  the  initial  SMD  close  to  the  nozzle  at  the  same  condition,  so 
that  the  change  in  SMD  with  distance  is  clearer.  Comparing  the  emulsified  versus 
neat  fuels  in  Figures  32-36,  the  trends  for  SMD's  are  similar  for  pressures  of  345  kPa 


■  Emulsified  Hexadecane 


FIGURE  31b.  EFFECT  OF  AIR  TEMPERATURE  ON  SMD  OF  EMULSIFIED  AND 
NEAT  HEXADECANE  AT  26  mm  FROM  NOZZLE,  AIR  PRESSURE  OF  448  kPa  (4.42  ATM) 
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FIGURE  31a.  COMPARISON  OF  RELATIVE  SMD  FOR  EMULSIFIED  AND  NEAT  HEXADECANE  AT 
CONSTANT  PRESSURE  (4.42  atm,  448  kPa)  AND  VARIOUS  TEMPERATURES 
(DELAVAN  WDA  3.8  LITERS/HR  45°  HOLLOW  CONE,  AP=483kPa).  SMDinitial  =  28 


FIGURE  32.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

HEXADECANE  644  K,  165  kPa 


HEXADECANE  644  K,  241  kPa 
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FIGURE  34.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

HEXADECANE  644  K,  345  kPa 


DISTANCE  FROM  NOZZLE  (MM) 

FIGURE  35.  CHANGE  OF  SMD  AND  N  WITH  DISTANCE  FROM  THE  NOZZLE 

HEXADECANE  644  K,  448  kPa 


and  below,  but  diverge  at  448  kPa  with  the  emulsified  fuels  showing  smaller  drop  sizes. 
These  air  conditions  allow  certain  steady-state  drop  temperatures  to  be  reached,  as 
shown  in  Table  4.  The  divergence  in  drop-size  distributions  between  345  kPa  and  448 
kPa  (at  a  fixed  air  temperature  of  644  K)  corresponds  to  drop  temperatures  in  the 
range  of  512  K  to  520  K.  Thus,  for  variations  in  air  pressure,  the  divergence  in  drop- 
size  distributions  for  emulsified  and  neat  fuels  occurs  in  the  range  of  drop  tempera¬ 
tures  from  roughly  512  K  to  520  K,  similar  to  the  range  of  503  K  to  520  K  reported 
above  for  the  divergence  with  varying  air  temperature.  Both  changes  in  air 
temperature  and  air  pressure  lead  to  a  divergence  in  drop-size  distributions  for 
emulsified  and  neat  hexadecane  at  drop  temperatures  of  approximately  515  K.  This  is 
a  reasonable  drop  temperature  to  expect  microexplosions  in  a  practical  system 
compared  with  the  upper  limit  of  539  K  in  a  pure  system. 

In  addition  to  the  comparative  results  for  emulsions  versus  neat  fuels,  the  drop- 
size  data  for  the  neat  fuels  alone  provide  some  interesting  and  unique  data  with  which 
to  compare  predictions  from  computer  models  for  spray  evaporation  at  elevated 
temperatures  and  pressures.  Such  comparisons  are  presented  in  a  later  section. 

Data  for  both  jet-A  and  hexadecane  indicate  that  emulsified  and  neat  fuels 
initially  atomize  to  approximately  the  same  drop  sizes  for  a  given  condition,  and 
evaporate  in  such  a  way  so  as  to  maintain  similar  drop  sizes  at  lower  air  temperatures 
and  pressures.  As  air  temperatures  and  pressures  increase,  the  drop  sizes  diverge  and 
the  emulsified  drops  are  consistently  smaller  downstream  of  the  nozzle.  For 
hexadecane,  this  divergence  occurs  at  a  predicted  drop  temperature  on  the  order  of 
515  K.  Although  these  results  are  consistent  with  what  would  be  expected  if 
microexplosions  occurred  for  the  emulsified  fuels,  it  is  again  possible  that  the  higher 
specific  heat  and  heat  of  vaporization  of  water  could  delay  the  evaporation  process  for 
emulsions  and  produce  similar  results. 

3.2  Two-Phase  Probe  Measurements 


Although  microexplosions  of  emulsified  drops  may  produce  similar  drop-size 
results  to  those  caused  by  reduced  evaporation  rates  of  the  small  drops  due  to  the 
heat-sink  effect  of  water,  these  two  phenomena  have  opposite  effects  on  the  spray 
evaporation  rates.  The  heat-sink  effect  would  reduce  evaporation  rates  while 
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microexplosions  would  enhance  the  evaporation  due  to  the  smaller  drop-sizes  pro¬ 
duced.  Therefore  a  measurement  of  the  comparative  evaporation  rates  of  the 
emulsified  and  neat  fuels  would  determine  the  relative  importances  of  the  micro¬ 
explosions  versus  the  heat  sink  effect. 

Such  measurements  of  the  fuel  vapor  concentration  across  the  spray  were 
performed  using  the  two-phase  probe  described  earlier.  Measurements  were  per¬ 
formed  at  two  distances  from  the  nozzle,  19  mm  and  57  mm,  at  temperatures  from 
atmospheric  to  650  K,  and  at  a  constant  pressure  of  448  kPa  (4.42  atm).  Profiles  for 
emulsified  and  neat  hexadecane  fuel  vapor  concentration  (hydrocarbon  concentration 
only)  at  57  mm  from  the  nozzle  are  shown  in  Figure  37  for  509  K  and  Figure  38  for 
646  K.  Similar  profiles  were  recorded  at  other  temperatures  at  57  mm,  and  at  19  mm 
from  the  nozzle.  Equal  volume  flowrates  were  used  for  the  emulsified  and  neat  fuels, 
so  that  the  hydrocarbon  concentration  for  complete  evaporation  of  the  emulsified  fuel 
would  be  only  about  80  percent  of  that  for  the  neat  fuel. 

The  profiles  for  57  mm  such  as  those  shown  in  Figures  37  and  38,  as  well  as  the 
profiles  at  19  mm,  were  integrated  across  the  whole  spray,  accounting  for  the  area 
weighting  factors  and  correcting  the  emulsified  fuels  for  the  20  percent  hydrocarbon 
displacement  by  water,  to  get  a  single  integrated  value  for  each  fuel  at  a  given 
condition.  These  integrated  hydrocarbon  concentrations  are  plotted  as  a  function  of 
air  temperature  in  Figure  39  for  the  19  mm  location  and  Figure  40  for  the  57  mm 
location.  (Actually,  the  effective  sampling  location  for  the  probe  is  somewhat 
upstream  of  the  physical  location  of  the  probe  tip  due  to  the  acceleration  of  the  flow 
into  the  probe,  but  no  estimate  is  included  to  account  for  this.) 

As  shown  in  Figures  39  and  40,  very  little,  if  any,  fuel  vapor  could  be  measured 
at  306  K,  in  agreement  with  the  low  vapor  pressure  of  hexadecane  at  this  temperature. 
However,  this  is  very  significant  in  that  it  indicates  that  the  aerodynamic  separation 
of  drops  from  the  vapor  sample  works  very  effectively,  even  in  these  finely  atomized 
sprays.  The  probe  itself  is  heated,  so  if  drops  were  present  in  the  vapor  sample,  they 
would  be  vaporized  and  the  measured  values  would  be  greater.  This  room  temperature 
(approximately)  measurement  also  indicates  that  the  heated  probe  is  not  affecting  the 
sample  by  vaporizing  the  spray  either  before  or  after  entering  the  probe  and  prior  to 
the  "vapor  only"  separation  point  25  mm  downstream  of  the  probe  tip.  The  data  in 
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Emulsified  Hexadecane 


FIGURE  39.  COMPARISON  OF  HYDROCARBON  VAPOR  CONCENTRATION  FOR  EMULSIFIED 

AND  NEAT  HEXADECANE  AT  19  mm 


FIGURE  *0.  COMPARISON  OF  HYDROCARBON  VAPOR  CONCENTRATION  FOR  EMULSIFIED 

AND  NEAT  HEXADECANE  AT  57  mm 


Figure  39  indicate  that  at  19  mm,  the  vaporization  of  the  neat  and  emulsified  fuels  is 
similar  within  the  precision  of  the  measurement.  At  57  mm,  the  emulsified  fuel  is 
vaporizing  faster  at  elevated  temperatures.  This  observation  lends  strong  support  to 
the  idea  that  sufficient  microexplosions  are  occurring  to  enhance  the  vaporization  rate 
beyond  any  reduction  caused  by  the  heat-sink  effect  of  the  water. 

3.3  Comparison  of  Computer  Predictions  and  Measured  Spray  Parameters 


The  data  for  the  evaporation  of  neat  hexadecane  at  elevated  temperatures  and 
pressures  provides  not  only  a  reference  for  the  emulsified  fuel  data,  but  also  a  set  of 
data  for  comparison  with  predictions  from  the  computer  model.  Diagnostics  have  been 
described  which  provided  a  measure  of  both  drop-size  distributions  as  a  function  of 
distance  from  the  nozzle  and  hydrocarbon  vapor  concentration.  The  computer  model 
described  earlier  in  this  report  was  used  to  predict  these  values  in  a  form  comparable 
with  the  experimental  results.  Results  for  predicted  versus  measured  changes  in  drop- 
size  distribution  are  presented  first,  followed  by  comparisons  for  the  fuel  vapor 
concentration. 

The  computer  model  does  not  predict  initial  drop-size  distribution  upon  atomiza¬ 
tion,  but  requires  this  as  an  input.  The  experimental  values  measured  at  16.4  mm  were 
used  as  a  calibration  point  and  different  parameters  for  the  Rosin-Ram mler  size- 
distributions  at  the  nozzle  exit  were  used  until  a  reasonable  match  was  obtained  at 
16.4  mm.  The  computet  model  and  experimental  values  were  not  matched  closer  to 
the  nozzle  because  of  the  sampling  bias  effects  of  the  laser  beam  discussed  previously. 
That  is,  too  close  to  the  nozzle,  the  large  drops  on  the  outer  periphery  cause  a 
significant  increase  in  the  measured  SMD,  an  effect  not  included  in  the  computer 
model.  Thus,  the  comparison  between  predicted  and  measured  size  distributions  should 
be  made  for  the  trends  of  SMD  and  N  as  a  function  of  distance  from  the  nozzle  and  air 
temperature. 


Predictions  from  the  computer  model  are  compared  with  experimental  results 
for  a  series  of  tests  at  44S  kPa  (4.42  atm),  and  a  temperature  of  450  K  in  Figure  41, 
509  K  in  Figure  42,  547  K  in  Figure  43,  593  K  in  Figure  44,  and  646  K  in  Figure  45. 
The  lines  in  each  figure  are  the  model  predictions  and  the  symbols  are  the 
experimental  values.  In  each  figure,  the  upper  line  and  data  are  for  the  Rosin- 
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FIGURE  44.  COMPARISON  OF  PREDICTED  AND  OBSERVED  SMD  AND  N  PARAMETER 
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FIGURE  45.  COMPARISON  OF  PREDICTED  AND  OBSERVED  SMD  AND  N  PARAMETER 

646  K,  448  kPa 


Rammier  N  parameter,  and  for  each  case,  the  predicted  and  measured  trends  ; 
fairly  well,  with  both  indicating  a  narrower  distribution  (larger  N)  going  downst 
from  the  nozzle,  corresponding  to  the  evaporation  and  removal  from  the  distrib 
of  the  smaller  drops.  Figure  45  for  646  K  does  show  a  discrepancy  between  pred 
and  measured  N  values,  but  the  measured  value  is  somewhat  in  error  due  to  I 
steering  problems.  At  the  highest  temperatures,  the  correction  technique  for  e 
due  to  steering  of  the  laser  beam  by  thermal  gradients  cannot  completely  correc 
light-scattering  signal,  with  the  results  that  a  systematic  error  is  introduced  ca 
the  measured  N  to  be  smaller  than  its  true  value  and  the  SMD  to  be  larger  tha 
true  value.  This  limits  the  application  of  this  instrument  at  some  point  downstrea 
the  nozzle. 

The  predicted  and  measured  SMD's  are  shown  in  the  lower  parts  of  Figures  4 
Notice  that  both  the  lines  for  the  computer  model  predictions  and  the  syn 
representing  the  measured  values  increase  more  rapidly  going  downstream  at 
higher  temperatures.  These  trends  seem  to  be  predicted  fairly  accurately,  with  i 
some  exception  in  Figure  45.  The  errors  in  SMD's  due  to  beam  steering  are  usually 
severe  than  the  errors  in  N's.  The  predicted  and  measured  SMD’s  . -;o  diverge  clo 
the  nozzle  due  to  the  sampling  effects  discussed  earlier,  which  also  cause  some  o 
differences  between  the  predicted  and  measured  values  at  the  lower  temperature 
particular,  the  measured  values  being  higher  than  predicted  close  to  the  nozzle 
lower  downstream. 

The  computer  model  was  also  used  to  predict  the  fuel  vapor  concentration  ii 
spray  cone,  and  these  values  were  used  for  comparison  with  the  two-phase  p 
measurements  for  the  hexadecane  spray.  The  comparison  at  19  mm  downstream  o 
nozzle  is  shown  in  Figure  46,  and  at  57  mm  downstream  in  Figure  47.  The  comp 
model  assumes  a  uniform  concentration  in  the  radial  direction  to  the  edge  of  the  s 
cone  as  shown  in  Figure  5.  It  is  encouraging  that  the  values  are  of  similar  magniti 
but  there  are  some  considerable  differences  that  are  due  to  both  simplT 
assumptions  in  the  model  and  questions  about  the  two-phase  probe.  These  re 
could  be  improved  by  a  measurement  of  the  liquid-phase  luel  concentration  along 
the  vapor  phase.  Such  a  measurement  is  possible  with  the  two-phase  probe,  and  \» 
allow  the  fraction  of  fuel  which  is  vaporized  to  be  calculated  more  accurately 
measurement  of  the  liquid  fuel  concentration  as  a  function  of  radial  displacer 
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FIGURE  46.  COMPARISON  OF  PREDICTED  AND  MEASURED  HYDROCARBON 
VAPOR  CONCENTRATION  AT  19  mm 
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FIGURE  47.  COMPARISON  OF  PREDICTED  AND  MEASURED  HYDROCARBON 
VAPOR  CONCENTRATION  AT  57  mm 


from  the  nozzle  axis  would  also  allow  a  check  of  the  trajectory  portion  of  the 
computer  model. 

The  detailed  comparison  presented  here  between  predicted  and  measured  drop- 
size  distributions  and  fuel-vapor  concentrations  for  fuel  sprays  into  elevated  tempera¬ 
ture  and  pressure  air  shows  the  promising  capabilities  of  these  diagnostic  techniques 
for  model  verification.  Verification  data  for  fuel  spray  models  is  severely  lacking,  and 
these  diagnostics  demonstrate  an  approach  to  model  verification. 

3.4  Drop-Size  Measurements  in  Burning  Sprays 

The  measurements  and  modeling  discussed  earlier  have  been  used  for  non-burning 
sprays  injected  into  uniform  temperature  air  flowing  coaxially  with  the  fuel  spray.  It 
was  also  of  interest  to  attempt  drop-size  measurements  in  a  burning  spray.  Because  of 
the  sharp  temperature  gradients  and  unsteadiness  of  the  combustion  process,  the 
results  would  be  difficult  to  interpret,  but  it  could  be  a  useful  experimental  diagnostic 
instrument.  Two  types  of  problems  are  present  when  making  drop-size  measurements 
in  burning  sprays.  First,  the  flame  produces  high  levels  of  visible  and  near-infrared 
radiation  coinciding  with  the  response  range  of  the  detectors  in  the  sizing  instrument, 
and  orders  of  magnitude  higher  in  intensity  than  the  desired  scattered  radiation  from 
the  HeNe  laser.  Secondly,  the  strong  gradients  in  temperature  and  air/fuel  vapor 
concentration  result  in  beam  steering  problems  even  worse  than  those  in  the  non¬ 
burning  evaporating  sprays.  Oust  as  for  sprays  evaporating  in  elevated-temperature 
air,  the  beam  steering  results  in  extraneous  amounts  of  light  being  recorded  on  the 
inner  detector  rings,  resulting  in  significant  errors  in  the  computed  drop-size  distribu¬ 
tion. 

Drop-size  measurements  were  attempted  with  JP5  in  emulsified  (20-percent 
water)  and  neat  form  in  two  different  types  of  flames.  First  the  disc-in-duct 
combustor  described  earlier  was  used.  Then  the  flame  was  stabilized  downstream  of 
the  nozzle  using  two  rectangular  flame  holders,  placed  near  the  top  and  bottom  of  the 
spray,  42  mm  downstream  of  the  nozzle  and  separated  by  25  mm. 

For  both  flame  configurations,  the  phase-sensitive  detection,  optical  filtering, 
and  restriction  of  the  field-of-view  of  the  detector  described  in  the  EXPERIMENTAL 


APPARATUS  section  were  effective  in  preventing  interference  from  the  flame 
radiation  in  the  scattered-light  measurement.  However,  the  extreme  thermal  gradi¬ 
ents  in  the  flame  regions  prevented  reasonable  data  from  being  obtained  due  to 
excessive  beam-steering  problems.  For  the  case  where  the  flame  was  stabilized 
downstream  of  the  spray  nozzle,  drop-size  data  were  obtained  from  the  nozzle  up  to 
the  beginning  of  the  visible  flame.  For  these  fairly  fine  and  dilute  sprays,  the 
evaporation  at  the  flame  front  is  so  rapid  and  the  temperature  gradients  so  steep  that 
no  meaningful  data  were  obtained  within  the  burning  portion  of  the  spray.  When  the 
flame  was  stabilized  downstream  of  the  nozzle,  the  measured  SMD  was  fairly  constant 
to  within  about  4  mm  of  the  visible  flame,  and  then  increased  from  there  to  the  edge 
of  the  flame  front  where  data  could  no  longer  be  acquired. 

It  appears  that  drop-size  measurements  in  burning  sprays  are  not  possible  with 
the  diffraction-based  instruments  unless  other  approaches  are  used  to  alleviate  the 
beam  steering  problems.  In  addition,  the  vaporization  is  so  rapid  that  very  good 
spatial  resolution  would  also  be  required  for  the  conditions  of  this  experiment. 
Techniques  to  overcome  the  other  problem  of  interference  of  flame  radiation  have 
been  developed  and  described  here. 


IV.  SUMMARY  AND  DISCUSSION 


Significant  progress  has  been  made  during  this  program  in  three  general  areas. 
First,  diagnostics  have  been  developed  for  measuring  spray  characteristics,  including 
drop-size  distributions  and  vaporized  fuel  concentrations,  in  elevated-temperature  and 
pressure-evaporating  sprays.  Second,  these  diagnostics  have  been  used  to  determine 
the  relative  atomization  and  evaporation  characteristics  of  fuels  emulsified  with 
water  and  their  neat  counterparts  in  polydisperse  high-velocity  fuel  sprays.  The 
purpose  of  these  comparative  measurements  has  been  to  determine  if  microexplosions 
can  occur  in  emulsified  fuels  with  sufficient  frequency  and  intensity  to  enhance  the 
overall  evaporation  rate  when  compared  with  neat  fuels.  Third,  these  experimental 
measurements,  when  combined  with  predicted  results  from  a  spray  trajectory  and 
evaporation  computer  model,  have  led  to  a  better  understanding  of  fuel  spray 
evaporation. 

The  data  presented  here  for  emulsified  and  neat  hexadecane  indicate  a  diver¬ 
gence  in  drop  size  distributions  at  a  predicted  drop  temperature  of  about  515  K.  The 
predicted  spontaneous  nucleation  temperature  of  water  in  hexadecane  is  about  539  K. 
However,  the  presence  of  particles,  bubbles,  or  dissolved  gases  would  lower  the  actual 
nucleation  temperature  below  the  predicted  value.  The  emulsions  used  in  these  tests 
were  ultrasonicated  to  keep  them  uniformly  mixed  before  entering  the  test  section, 
and  bubbles  may  have  been  created  within  the  water  or  at  the  water/surfactant/fuel 
interface.  Thus,  the  observed  divergence  in  the  drop-size  distributions  of  neat  and 
emulsified  fuels  occurs  at  a  predicted  drop  temperature  reasonably  close  to  the 
expected  nucleation  temperature  of  water. 

The  trends  for  drop-size  distributions  observed  for  emulsified  and  neat  jet-A 
were  qualitatively  similar  to  those  seen  for  the  hexadecane  blends.  Because  of  the 
multicomponent  nature  of  jet-A,  drop  temperatures  are  more  difficult  to  estimate  and 
no  attempt  was  made  to  calculate  those  temperatures  here. 

The  measurements  with  the  two-phase  probe  allow  a  distinction  between  drop- 
size  differences  caused  by  the  heat-sink  effect  of  water  (due  to  its  large  sensible  and 
latent  heat  of  vaporization)  and  microexplosions.  The  heat-sink  effect  will  reduce  the 
evaporation  rate  of  emulsions  relative  to  the  neat  fuel,  while  the  microexplosions 


should  increase  it.  Vapor-phase  measurements  show  that  the  emulsified  and  neat  fuels 
initially  evaporate  at  about  the  same  rate,  but  further  from  the  nozzle  and  at  higher 
temperatures  and  pressures  the  emulsified  fuel  evaporates  about  40  percent  faster. 
Microexplosions  are  the  probable  explanation  for  these  experimental  results.  Thus, 
even  in  a  turbulent  high-pressure  air  stream  where  internal  drop  mixing  is  enhanced, 
rrucroexplosions  can  occur,  creating  smaller  droplets  to  enhance  the  overall  evapora¬ 
tion  rate;  the  effect  is  apparently  sufficient  to  overcome  the  reduction  in  evaporation 
rate  due  to  the  high  heat  capacity  of  water.  The  effect  of  air  pressure  on 
microexplosions  was  discussed  earlier,  but  it  should  also  be  mentioned  that  while 
increasing  pressure  enhances  the  probability  of  microexplosions,  the  intensity  of  the 
disruption  is  expected  to  be  reduced .(23) 

This  research  does  not  resolve  the  question  about  why  emulsification  of  fuels 
with  water  reduces  soot  formation,  but  it  does  verify  that  microexplosions  cannot  be 
excluded  as  an  explanation  even  in  turbulent  air  streams. 

Aside  from  the  questions  about  emulsified  fuels,  the  diagnostics  and  model 
discussed  here  have  application  to  spray  evaporation  problems  in  general.  Faeth's^D 
review  article  on  spray  models  indicates  that  there  is  a  scarcity  of  experimental 
confirmation  of  evaporation  codes.  In  order  to  verify  spray  models,  several  diagnostic 
devices  are  required,  and  two  have  been  described  here. 

Measurements  with  the  drop-sizing  instrument  have  shown  that  nozzle  cali¬ 
brations  obtained  at  atmospheric  conditions  cannot  be  assumed  to  be  an  accurate 
estimate  of  nozzle  performance  in  a  combustor.  The  drop-sizing  instrument  as 
modified  for  this  work  can  be  used  to  obtain  initial  spray  sizes  of  many  types  of  sprays 
at  elevated  temperatures  and  pressures.  However,  laser  beam  steering  caused  by 
temperature  and  concentration  gradients  limits  its  application  at  some  point  down¬ 
stream  of  the  nozzle.  For  the  experiments  described  here,  data  could  be  obtained  only 
for  the  first  20  percent  or  so  of  the  evaporation  process. 

Two-phase  or  phase-discriminating  probes  offer  the  capability  of  determining 
fuel  vaporization  rates  within  a  drop-laden  fuel  spray.  The  aerodynamics  separation 
appears  to  work  well  even  for  the  fine  sprays  (SMD  =  30  micrometers)  used  in  this 
work,  as  indicated  by  the  room  temperature  tests  with  hexadecane.  The  probe  did 


seem  to  exhibit  some  hysterisis  with  changes  in  air  temperatures  that  is  not  yet 
understood.  It  would  be  desirable  to  test  and  confirm  the  operation  of  a  two-phase 
probe  using  an  independent  technique  such  as  a  non-intrusive  optical  absorption 
method. 
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V.  CONCLUSIONS 


Drop-size  distributions  and  fuel  vapor  concentrations  have  been  determined  at  a 
variety  of  elevated  temperatures  and  pressures  within  sprays  of  jet-A  and  hexadecane 
and  their  emulsions  with  20  percent  water.  These  results  indicate  that  microexplo¬ 
sions  occur  in  the  emulsified  fuel  at  some  conditions  in  spite  of  the  high  drop-to-air 
Reynolds  numbers  and  the  air  turbulence  which  serve  to  enhance  internal  drop 
circulation  and  reduce  the  probability  of  microexplosions. 

Modifications  have  been  developed  and  described  which,  for  the  first  time,  allow 
the  Malvern  laser-diffraction  drop-sizing  instrument  to  be  used  in  high-temperature 
rapidly-evaporating  sprays.  Modifications  are  also  described  for  this  instrument  to 
overcome  problems  associated  with  high  background  light  levels  such  as  flame 
radiation.  However,  in  the  case  of  flames,  the  beam-steering  problems  become  very 
severe.  A  calibration  procedure  for  the  Malvern  has  been  developed  and  published 
which  results  in  agreement  with  calibration  verification  reticles  to  within  the 
accuracy  of  the  reticles.  This  procedure,  described  in  Reference  3  of  Appendix  B, 
should  alleviate  the  differences  between  instruments  which  have  been  so  troublesome 
in  the  past. 

Results  have  also  been  reported  for  a  two-phase  probe.  This  appears  to  be  a 
powerful  diagnostic  instrument  for  spray  model  verification,  but  more  work  needs  to 
be  done  to  resolve  questions  concerning  its  operation. 

The  computer  model  for  spray  trajectory  and  evaporation  described  here  predicts 
many  of  the  trends  observed  experimentally  for  both  the  drop-size  data  and  fuel  vapor 
concentration  data. 

The  air  pressure  and  temperature  have  a  considerable  effect  upon  initial 
atomization  of  the  spray,  with  increasing  air  density  improving  atomization.  This  has 
important  consequences  when  using  drop-size  calibrations  obtained  at  atmospheric 
conditions  to  predict  performance  of  gas-turbine  combustors  at  elevated  pressure,  or 
at  the  low  air  densities  typical  of  altitude-relight  conditions. 


VI.  RECOMMENDATIONS 


The  results  from  this  program  suggest  several  other  interesting  areas  for 

applications  of  techniques  developed  here. 

1.  The  diagnostics  which  were  developed  in  this  program  to  study  secondary 
atomization  in  emulsified  fuel  could  be  used  to  examine  secondary  atomization 
in  azide/fuel  mixtures,  or  other  blends  designed  for  enhanced  atomization. 

2.  The  comparisons  of  predictions  from  the  spray  model  with  measured  quantities 
show  some  interesting  capabilities  for  the  model.  The  computer  model  should  be 
expanded  to  include  more  of  the  sampling  effects  and  the  spray/ air  momentum 
exchange,  and  a  more  thorough  comparison  made  with  existing  experimental  data 
as  well  as  newly  acquired  data. 

3.  The  two-phase  probe  tests  were  limited  in  scope,  but  indicated  a  powerful 
technique  for  spray  model  verification.  This  technique  should  be  expanded  to 
include  liquid-fuel  analysis.  Also,  the  results  should  be  determined  independently 
by  an  optical  absorption  technique  to  develop  more  confidence  in  the  probe. 

4.  It  is  shown  in  this  report  that  increasing  air  density  improves  atomization.  The 
effect  of  density  and  temperature  need  to  be  documented  more  generally  for 
subatmospheric  and  elevated  pressures.  This  has  important  implications  for 
modeling  gas-turbine  combustion  and  altitude-relight  problems. 

5.  Dense  fuel  sprays  cause  errors  for  the  drop-sizing  instrument  due  to  multiple¬ 
scattering  of  photons  off  more  than  one  drop.  This  effect  has  been  examined 
experimentally  and  some  empirical  correction  factors  developed,  but  a  theoreti¬ 
cal  study  is  desirable  to  develop  a  more  general  correction  scheme. 
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APreNDIX  A 


FUEL  SPRAY  COMPUTER  MODEL  -  MATHEMATICS 


For  purposes  of  discussion,  the  model  is  divided  into  three  parts 
called  (1)  Thermodynamics,  (2)  Aerodynamics,  and  (3)  Drop-Size 
Distribution.  Each  section  is  discussed  in  order.  An  overall  flow  chart 
is  shown  in  Figure  A-l . 

Thermodynamics 

This  part  of  the  model  computes  the  heat-up  of  the  drop,  the  final 
steady-state  temperature,  and  the  fuel  and  air  properties  necessary  for 
those  calculations.  This  part  of  the  model  is  identical  to  that 
described  by  Chin  and  Lefebvre  (Ref.  Al  and  A2)  except  that  some  of  the 
fuel  property  data  were  taken  from  the  American  Petroleum  Institute  Data 
Book  (Ref.  A3).  The  results  from  this  model  are  very  similar  to  those 
presented  by  Chin  and  Lefebvre  when  the  model  is  used  for  quiescent  fuel 
sprays.  The  Integration  of  this  part  of  the  model  with  the  other 
sections  to  be  described  allows  the  examination  of  sprays  in  situations 
where  both  spray  and  ai.  are  moving  with  a  nonzero  relative  velocity. 
The  mathematics  are  outlined  below,  but  a  more  detailed  derivation  is 
given  by  Chin  and  Lefebvre  (Ref.  AL  and  A2)  wh'ch  in  turn  Is  based  on  the 
theories  described  by  Spalding  (Ref.  A4 ; . 

The  steady-state  temperature  Is  determined  as  follows.  The  mass 
transfer  rate,  m,  from  a  drop  Is  given  by, 

mp  2nL)  (kg/c^g)  In  (1  +  B)  (Al) 

where  D  is  the  diameter,  kg  Is  the  thermal  conductivity  of  the  gas,  cp,g 
is  the  specific  heat  at  constant  pressure,  and  B  is  either  Bvj,  the  mass 
diffusion  transfer  number  If  mass  diffusion  is  controlling,  or  B-p,  the 
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MB 


Input: 

Fuel  properties 
Air  properties 
Nozzle  characteristics 
Drop  size  distribution 


_ ^ 

Assume  steady  -  state 

.  Iteration 

Update  fuel  and  air 

drop  temperature  and 

< - > 

properties  for  assumed 

check  =  B-p  ? 

steady- state  temperature 

Final  steady  -  state 
drop  temperature  when 
BM  ~  Bj 


Calculate  drop  temperature  during 
heat-up  phase  accounting  for  changing 
fuel  and  air  properties 


If  heat-up  95%  complete  use 
steady-state  temperature 


Calculate  mass  loss  rate 
and  new  drop  size 


Calculate  drag,  acceleration, 
velocity  and  position 


Check  to  see  if  drop  has 
arrived  at  desired  position 


»V  :  m 


N0  I  Change  drop  size 
All  drops  calc.? 


When  all  drops  have  arrived  at  a 
desired  position,  calculate  drop 
size  distribution 


Update  desired  position 


thermal  diffusion  transfer  number  if  heat  transfer  is  controlling  and  the 
Lewis  number  is  assumed  to  be  unity.  is  defined  as 

BM  =  (YF,s  -  YF>0o)  /  (1  -  YF,s)  (A2) 

where  Ypjg  and  Yp^  are  the  fuel  mass  fractions  at  the  drop  surface  and 
in  the  ambient  air,  respectively,  and  Yp>g  is 

YF,s  =  Pf.s  %  /  (PF.s  MF  +  (P  -  PFfs)  MA)  (A3) 

where  Pp>8  Is  the  fuel  vapor  pressure  at  the  drop  surface,  P  is  the 
ambient  pressure,  and  Mp  and  MA  are  the  molecular  weights  of  fuel  and 
air,  respectively.  The  vapor  pressure  was  calculated  using  the  method 
recommended  by  the  American  Petroleum  Institute  Data  Book  (Ref.  A3) 
called  Procedure  5A1.10,  "Vapor  Pressures  of  Pure  Hydrocarbons." 

The  thermal  diffusion  transfer  number  B-p  for  an  evaporating  drop  is 

bT  -  cP)g  (Tm  -  Ts)  /  L  (A4) 

where  Too  and  Ts  are  the  ambient  air  temperature  and  drop  surface 
temperature  and  L  is  the  latent  heat  of  fuel  vaporization  corrected  from 
the  normal  boiling  point  to  the  actual  surface  temperature. 

The  accuracy  of  Equation  (AJ_)  is  very  dependent  on  the  choice  of 
values  for  kg  and  Cp^.  As  recommended  by  Chin  and  Lefebvre  (Ref.  Al)  a 
reference  temperature  was  chosen  as  the  drop  surface  temperature  plus  1/3 
of  the  difference  between  ambient  air  and  surface  temperatures. 
Similarly  a  reference  value  of  the  fuel  vapor  mass  fraction  was  taken  as 
the  value  just  outside  the  surface  plus  1/3  of  the  difference  between  the 
value  at  infinity  and  just  outside  the  surface.  Using  these  reference 
conditions  denoted  by  ”r"  the  specific  heat  of  the  gas  is  given  by, 


Cp.g  =  YA,r  (cp,A  at  Tr)  +  Yp>r  (cp,v  at  Tr)  (A5) 

and  the  thermal  conductivity  by, 

kg  •  YA,r  (kA  at  Tr)  +  YF,r  (kV  at  Tr)  (A6) 


At  steady-state  conditions,  ■  Bf  or  from  Equations  (A2)  and  (A4) 
yF,s  -  Ypj0O  _  Cp,g  (T®  -  Ts>st) 


(AT) 


1-VF.s  L 

where  Ts>at  is  the  desired  surface  temperature  under  steady-state 
conditions.  Assuming  that  the  fuel  vapor  concentration  in  the  ambient 
air  is  zero  (Ypj00  ■  0)  and  substituting  for  Yp>a  from  Equation  (A3^)  then 
(A7)  becomes 

P  Mp  L 

pF,s 


mA  cp,g  (T<=°  *  Ts,st) 


-1  =  0 


(AS) 


In  solving  Equation  (A8)  to  determine  the  steady-state  surface 
temperature  TajSt,  some  of  the  variables  can  be  specified,  P,  Mp,  MA,  and 
To,  but  the  remaining  three  variables  are  functions  of  temperature,  with 
L  being  related  to  the  heat  of  vaporization  at  the  normal  boiling  point 


^T.bn 


k  =  kj.bn  ["  (Tcr  *  Ts)  /  (Tcr  -  Tbn)  J 


0.38 


(A9) 


where  Tcr  is  the  critical  temperature  and  Tbn  Is  the  normal  boiling  point 
of  the  fuel.  In  order  to  solve  (A8)  it  is  necessary  to  first  assume  a 
value  for  Ta>at.  That  temperature  is  used  to  calculate  a  reference 
temperature  Tr  which  in  turn  is  used  to  calculate  those  quantities 
specified  in  Equations  and  (A6)  used  to  calculate  Cptg  and  kg,  and 
to  calculate  L  as  specified  in  (A9).  If  the  assumed  value  of  Ta>st  is 
too  low,  the  left  side  of  Equation  (A8)  will  be  positive.  An  iterative 
procedure  must  be  used  to  continue  specifying  TS|gt  until  Equation  (A8) 
Is  satisfied  within  the  desired  accuracy.  That  procedure  results  In  a 
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determination  of  the  steady  state  surface  temperature  Tg>gt  and  the 
transfer  number  The  steady  state  quiescent  evaporation  constant 
^q  st  can  then  be  determined  from. 


Aq.st  =  8  kg  In  (1+B)  /  *F  cp>g 


(A10) 


where  is  the  fuel  density.  Frossling  (Ref.  A5)  has  shown  that 

convective  effects  can  be  accounted  for  in  the  case  where  heat  transfer 


rates  are  controlling  by. 


AC,st  =  Xq,st  (1  +  0.276  ReD0,5  Prg0*33) 


(All) 


where  Xc,st  Is  the  steady-state  evaporation  constant  corrected  for 
convective  effects,  Prg  is  the  Prandtl  number  for  the  gas,  and  Rej)  is  the 
Reynolds  number  using  the  relative  velocity  between  the  drop  and  the  gas. 
This  velocity  should  include  a  fluctuating  component  in  the  case  of 
turbulent  flow.  Turbulence  intensities  were  not  measured  in  these 
experiments,  but  a  value  of  20  percent  of  the  axial  velocity  was  used  as 
an  estimate.  The  effective  velocity  used  to  calculate  Rep  for  the 
evaporation  calculation  was  taken  as  the  sura  of  the  fluctuating  component 
(20  percent  of  the  air  velocity)  plus  the  difference  between  the  drop  and 
air  velocities.  These  computations  specify  the  steady-state  properties 
of  the  drop. 


The  evaporation  which  occurs  during  the  drop  heat-up  period  is 
significant  in  many  practical  situations,  with  the  drop  being  completely 
evaporated  before  reaching  steady-state  conditions  is  some  instances. 
Chin  and  Lefebvre  (Ref.  A2)  have  predicted  that  high  air  pressures  and 
convective  effects,  which  are  often  absent  in  laboratory  experiments, 
both  tend  to  Increase  the  relative  Importance  of  the  heat-up  period 
relative  to  the  steady-state  phase  of  drop  evaporation.  Because  the 
experiments  conducted  in  this  study  were  of  real  sprays  at  high  pressures 
and  relative  velocities,  it  was  necessary  to  model  the  heat-up  period  in 
detail. 


Chin  and  Lefebvre  (Ref.  A2)  have  shown  that  the  rate  of  change  of 
drop  surface  temperature  is  given  by, 

dTs  /  dt  =  (ri»F  L  /  cp,F  m)  (  (BT  /  BM)  -  1  )  (A12) 

where  mp  is  specified  by  Equation  (AO  with  B-Bfl  during  the  heat-up 
period,  and  m  is  the  drop  mass. 


m  =  (x/6)  Pp  ^  (A13) 

The  change  of  drop  size  with  time  is  given  by  (Ref.  A2), 

dD  /  dt  =  -4  kg  In  (1+BM)  /  ^p  cp>g  D  (A14) 

The  drop  temperature  asymptotically  approaches  the  steady-state 
temperature,  and  it  was  assumed  that  when  the  temperature  had  risen  to  95 
percent  of  the  difference  between  the  initial  fuel  temperature  and  the 
steady-state  temperature,  that  the  heat-up  calculations  could  be 
terminated  and  steady-state  properties  used. 

Aerodynamics 

The  approach  used  was  a  much  more  simplified  version  of  the  one 
developed  at  the  University  of  Sheffield  and  described  in  several  papers, 
e.g.,  Boyson  and  Swithenbank  (Ref.  A6).  A  cylindrical  coordinate  system 
was  used  where  x  is  the  axial  distance,  y  the  radial  distance,  and  z  the 
angular  position,  with  corresponding  velocities  in  the  axial  (u),  radial 
(v),  and  tangential  (w)  directions.  For  this  work,  the  tangential 
velocities  were  assumed  zero. 

The  equations  of  motion  of  a  particle  (drop)  neglecting  all  forces 
except  drag,  are. 


-  F  (wp  -  w„) 


(A  17) 


*  r 


where  the  ”p"  subscripts  refer  to  the  particle  and  the  "<»”  subscript 
refers  to  the  free  stream,  and  F  is  given  by, 


F  = 


(18  Mg  /  Pp 


)  (CD  Re  /  24) 


(A18) 


where  Mg  Is  the  gas  viscosity,  Pp  Is  the  particle  density.  Dp  is  the 
particle  diameter,  Cp  Is  the  drag  coefficient,  and  Re  is  the  Reynolds 
number  defined  as, 


Re 


DP 


(A19) 


The  drag  coefficient,  Cp,  is  given  by  (Ref.  A7), 


Cd  =  27  Re"0-84 

0  Re  -  80 

(A20) 

CD  =  0.271  Re0-217 

80  <  Re  *  104 

(A21) 

CD  «  2 

104  >  Re 

(A22) 

The  equations  of  motion  (A15  to  A17 )  are  solved  numerically  using  a  step 
sire  of  l  to  10  ^s,  and  the  equations  of  trajectory  are  solved  In  a 

similar  manner. 
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These  equations  describe  the  trajectory  of  a  particle  in  a  gas  stream. 
For  each  Iteration  In  the  trajectory  calculation,  the  drop  size  and 
temperature  are  updated  using  the  procedure  described  in  the 
Thermodynamics  section. 

At  the  end  of  each  Iteration  through  the  trajectory  calculations, 
the  current  calculated  axial  position  of  the  particle  is  checked  against 
a  target  value.  When  the  particle  reaches  that  target  value  for  axial 
position,  the  position,  transit  time,  size,  temperature,  and  other 
parameters  for  that  drop  are  frozen  and  the  calculations  are  repeated  for 
the  next  larger  drop  size.  After  all  drops  reach  a  targeted  axial 
position,  corresponding  to  a  measurement  location,  the  drop  size 
distribution  Is  calculated,  important  data  are  printed  out,  the  target 
position  Is  moved  downstream,  and  the  calculations  are  repeated. 

Drop  Size  Distribution 

The  computer  model  makes  all  calculations  on  individual  drops  of 
certain  sizes,  but  the  Instrumentation  used  to  measure  the  spray  works 
best  in  characterizing  the  drop  sizes  by  two  parameters  which  specify  a 
size  and  width  of  the  distribution.  Three  distributions  may  be  selected: 
normal,  log-normal,  and  Rosin-Ramraler.  The  Ros  in-Rammler  distribution 
has  provided  the  best  fit  of  the  sprays  studied  and  has  been  used 
throughout  this  program.  In  order  to  compare  the  predictions  from  the 
computer  model  with  the  experimental  results,  it  is  necessary  to  convert 
the  computer  predicted  drop  sizes  at  a  given  location  into  an  equivalent 
set  of  two  parameters  specifying  the  Rosin-Raramler  distribution. 
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If  R  represents  the  fraction  of  the  liquid  being  sprayed  contained 
in  drops  larger  than  diameter  D,  then  for  the  Rosin-Rammler  distribution 
(Ref.  A8), 

R  =  exp  (-(D  /  i)N  )  (A26) 

where  x  represents  a  size,  and  N  specifies  the  width  of  the  distribution. 
It  is  convenient  to  specify  a  spray  by  a  single  parameter  representing  an 
"average"  drop  size.  However,  a  straight  numerical  average  is  heavily 
weighted  towards  the  smallest  drops  which  are  extremely  plentiful  but 
which  contain  a  very  small  fraction  of  the  total  volume  of  the  liquid. 
For  combustion  processes,  the  surf ace-to-volume  ratio  is  important  to 
evaporation,  so  an  average  drop  size  may  be  chosen  which  has  a  surface- 
to-volume  ratio  representative  of  the  actual  spray,  and  such  an  average 
is  called  the  surface-volume  mean  diameter  or  Sauter  mean  diameter  (SMD). 
The  SMD  is  related  to  the  Rosin-Rammler  parameters  by  (Ref.  A8), 

SMD  =  x  /  r  (1—1 /N)  N>1  (A“7) 

Specifying  x  and  N  allows  the  calculation  of  the  SMD ,  or  specifying  SMD 
and  N  determines  x. 

An  initial  set  of  drop  sizes  distributed  approximately  exponentially 
in  size  (i.e.,  evenly  spaced  when  plotted  as  In  D)  and  covering  the  range 
of  about  6(jra  to  560prr  (corresponding  to  the  response  range  of  the  drop 
sizing  instrument)  with  16  to  21  different  sizes  is  used  in  the  model. 
Assuming  size  classes  bounded  by  these  drops  as  end  points,  the  initial 
fraction  of  liquid  contained  in  drops  larger  than  a  certain  size  class  is 
given  by  Equation  ( A26 ) . 

Ri  =  exp  (-(D0>i  /  x)N)  l  <  i  <n  (A28) 

Ro  5  1 
Rn+1  =  0 


where  "n”  is  the  number  of  drops,  D0  is  the  initial  drop  diameter,  and 
the  initial  fraction  of  liquid  in  each  size  class  is. 


Fo,i  =  Ri  -  Ri+1 


Fo,-l  =0 


o  <  i  <  n 


(A29) 


After  evaporation  begins,  the  small  drops  begin  evaporating  quickly 
while  the  large  ones  evaporate  slowly,  changing  the  fraction  of  liquid  in 
the  different  size  classes  unevenly.  Denoting  the  smallest  nonevapo rated 
drop  size  by  "k”,  the  fraction  of  liquid  remaining  in  any  size  class  is, 


Fo,i  |(Di  /  Do,i)3  /  “p.i  +  C^i+i  /  D0>i+i)3  /  Up^  +  i  |/2 

n  7  "  '  '  ■  1  r 

51  F°>j  /  D°»j^3  /  up,i  +  (Dj  +  1  /  D0>j  +  i)3  /  Upj  +  1  |/2 

j=k 


(A  50) 


where 

D0  5  D1 

UP,0  2  upfl 

Dn+1  -  Dn 
up,n+l  =  “p.n 


k  <  i  <n 
(o<  k<n) 


and  upjj  is  the  axial  velocity  of  the  ith  drop,  is  the  instantaneous 
drop  size,  and  D0>i  is  the  initial  drop  size  of  the  ith  drop.  Each  size 
class  is  characterized  by  the  drops  making  up  the  end  points,  and  the 
initial  fraction  of  liquid  in  a  size  class  is  modified  by  the  average 
loss  of  volume  of  the  end  point  drops.  The  velocity  term  is  added  to 
account  for  the  fact  that  as  the  drops  slow  down  they  increase  their 
relative  concentration,  and  thus  their  weighting  factor,  in  the  sample 
vo lume . 
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The  cumulative  fraction  of  liquid  in  all  size  classes  including  and 
larger  than  the  ith  class  is. 


=£  Fi 


k  <  i  <n 


This  instantaneous  value  of  R^  is  related  to  the  drop  sizes  and 
Rosin-Rammler  parameters  by  Equation  (A28 )  with  the  initial  drop  size 
D0ji  replaced  by  the  instantaneous  value  D^, 


Ri  =  exp  (-(Di  /  x)‘  )  k  <  i  <  n 


(A32) 


Taking  the  natural  logarithm  of  both  sides  twice  and  excluding  i=k 
(Ri-1),  Equation  (A32)  becomes, 

lne( lne( 1  /  Ri))  =  N  lne  Di  -  N  lne  x  k+1  <  i  <  n  (A32 

This  has  the  form  of  the  equation  of  a  straight  line,  y=mx+b,  if  the 
following  definitions  are  used, 

Y  -  lne( lne( 1/R) ) 

m  ■  N 


x  *  lne  D 


b  -  -N  lne  x 

Thus,  by  determining  a  least  squares  fit  of  the  straight  line 
through  the  data  lne( lne( 1 /Rj ) )  versus  lneDi,  the  Rosin-Rammler 
parameters  are  given  by, 

N  =  m  (the  slope)  (A55 


x  =  exp  (-b  /  N) 


A  standard  routine  is  used  to  perform  the  least  squares  operation  and  the 
RoslnrRammler  parameters  and  the  SMD  (from  Equation  (A27 ) )  are  determined 
at  each  target  value  of  the  axial  location  corresponding  to  the  position 
where  experimental  data  are  obtained.  Although  the  Initial  distribution 
at  the  nozzle  Is  an  ideal  Ros ln-Rammle r  distribution,  the  distribution 
downstream  does  not  correspond  exactly  to  the  Ror  Rammler  equation  due 
to  the  different  evaporation  rates  for  the  difft  it  sized  drops,  and  the 
degree  of  fit  Is  determined  by  the  correlat  '  coefficient  of  the 
straight  line  through  the  computed  data. 
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Particle  Size  Measurement  Techniques,  and  accepted  for  publication  in  Liquid 
Particle  Size  Measurement  Techniques,  edited  by  J.M.  Tishkoff,  to  appear  1984. 

5.  "Atomization  and  Evaporation  of  Emulsified  and  Neat  Jet  Fuel  in  Real  Sprays," 
by  L.G.  Dodge  and  C.A.  Moses,  presented  at  the  Fall  1982  meeting  of  the 
Western  States  Section  of  the  Combustion  Institute. 


At  least  one  further  paper  describing  the  computer  spray  model  is  anticipated. 
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